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Abstract 



New data on the production of charged kaons in p+p interactions are presented. The data 
come from a sample of 4.8 million inelastic events obtained with the NA49 detector at 
the CERN SPS at 158 GeV/c beam momentum. The kaons are identified by energy loss 
in a large TPC tracking system. Inclusive invariant cross sections are obtained in inter- 
vals from to 1 .7 GeV/c in transverse momentum and from to 0.5 in Feynman x. Using 
these data as a reference, a new evaluation of the energy dependence of kaon production, 
including neutral kaons, is conducted over a range from 3 GeV to p+p colUder energies. 



1 Introduction 



Following the detailed investigation of inclusive pion lU and baryon |l2l production in 
p+p interactions, the present paper concentrates on the study of charged kaons. It thus completes 
a series of publications aimed at the exploration of final state hadrons in p+p collisions by using 
a new set of high precision data from the NA49 detector at the CERN SPS [|3|. The data have 
been obtained at a beam momentum of 158 GeV/c corresponding to a center-of-mass system 
(cms) energy of 17.2 GeV. This matches the highest momentum per nucleon obtainable with 
lead beams at the SPS, permitting the direct comparison of elementary and nuclear reactions. 
In addition, the chosen cms energy marks, concerning kaon production, the transition from 
threshold-dominated effects with strong s-dependences to the more gentle approach to higher 
energies where scaling concepts become worth investigating. On the other hand the character- 
istic differences between K+ and production which are directly related to the underlying 
production mechanisms, as for instance associate kaon-i-hyperon versus K-i-K pair production, 
are still well developed at SPS energy. They are manifest in the strong evolution of the K+/K~ 
ratio as a function of the kinematic variables. One of the aims of this paper is in addition the 
attempt to put the available results from other experiments into perspective with the present data 
in order to come to a quantitative evaluation of the experimental situation. 

A critical assessment of the complete s-dependence of kaon production seems the more 
indicated as its evolution in heavy ion interactions, especially in relation to pions, is promul- 
gated since about two decades as a signature of "new" physics by the creation of a deconfined 
state of matter in these interactions. As all claims of this nature have to rely completely on 
a comparison with elementary collisions, the detailed study of the behaviour of kaon produc- 
tion in p+p reactions from threshold up to RHIC and collider energies should be regarded as 
a necessity in particular as the last global evaluation of this type dates back by more than 30 
years [4J. A complete coverage of phase space, as far as a comparison of different experiments 
is concerned, is made possible in this paper, as compared to pions [IJ and baryons [2J, by the 
fact that there is no concern about feed-down corrections from weak hyperon decays, with the 
exception of Q decay which is negligible for all practical purposes. 

This paper is arranged in the same fashion as the preceding publications [|Tl|2|. A sum- 
mary of the phase space coverage of the available data from other experiments in Sect. |2] is 
followed by a short presentation of the NA49 experiment, its acceptance coverage and the cor- 
responding binning scheme in Sect. [3l Section |4] gives details on the particle identification via 
energy loss measurement as they are specific to the problem of kaon yield extraction. The eval- 
uation of the inclusive cross sections and of the necessary corrections is described in Sect. [51 
followed by the data presentation including a detailed data interpolation scheme in Sect. [6l 
K+/K^, K/vr and K/baryon ratios are presented in Sect.|7J A first step of data comparison with 
data in the SPS/Fermilab energy range is taken in Sect. [8l Section [9] deals with the data inte- 
grated over transverse momentum and the total measured kaon yields. The data comparison is 
extended, in a second step, over the range from y/s ~ 3 to ISR, RHIC and p+p collider energies 
in Sect.[lOl Section [TT] concentrates on an evaluation of K^. yields in relation to charged kaons 
and on a discussion of total kaon multiplicities as a function of y/s. A comment on the influ- 
ence of resonance decay on the observed patterns of pr and s dependence is given in Sect. [121 
In Sect. [131 a global overview of charged and neutral kaon yields as they result from the study 
of s-dependence in this paper is presented, both for the px integrated invariant yields at xp = 
and for the total kaon multiplicities. A summary of results and conclusion is given in Sect. [I4l 
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2 The experimental situation 



This paper considers the double differential inclusive cross sections of identified charged 



kaons, 



(1) 



as a function of the phase space variables defined as transverse momentum px and reduced 
longitudinal momentum 



Xf 



PL 

v/i/2 



(2) 



where pi denotes the longitudinal momentum component in the cms. 

If the phase space coverage of the existing data has been shown to be incomplete and 
partially incompatible for pion and baryon production in the preceding publications OlEl, the 
situation is even more unsatisfactory for charged kaons. A wide range of data covering essen- 
tially the complete energy range from kaon threshold via the PS and AGS up to the ISR and 
RHIC energy has been considered here. One advantage concerning the data comparison for 
kaons is the absence of feed-down from weak decays with the exception of Vt^ decay which can 
be safely neglected at least up to ISR energies. An overview of the available data sets is given 
in Fig. [T]for K+ and Fig. [2] for K~ in the x p/pr plane. 
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Figure 1: Phase space coverage of the existing K+ data: a) Cosmotron/PPA [I51|6l, b) PS/AGS 
BTHia, c) Serpukhov [[III, d) SPS/Fermilab [UMl, e) ISR [[IMS, f) RHIC [[23H301, g) NA49 



The sub-panels a) through g) show successively the energy ranges of the Cosmotron/PPA 
§M, PS/AGS [L-m, Serpukhov [11], SPS/Fermilab [UMl, ISR [15ji22] and RHIC [[23H301 
accelerators in comparison to the new data from NA49. The scarcity of data in the important 
intermediate energy range around ^/s ~ 10 GeV and the general lack of coverage in the low- 
Pt and low-xp regions are clearly visible. The coverage of the NA49 data, Figs.[T^ and[2k, is 
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Figure 2: Phase space coverage of the existing K data: a) Cosmotron/PPA, b) PS/ AGS [TTHTOI. 
c) Serpukhov LUJ, d) SPS/Fermilab [EMI, e) ISR [[l5Hni|2ffl2a, f) RHIC [EMOl, g) NA49 



essentially only limited by counting statistics towards high px and by limitations concerning 
particle identification towards high xp, in particular for K+, see Sect. |4]below. 

The task of establishing data consistency over the wide range of energies considered 
here is a particularly ardent one for kaons, as will be shown in the data comparison, see Sects. [8] 
and[10]below. This concerns especially any attempt at establishing total integrated yields where 
the existing efforts evidently suffer from a gross under-estimation of systematic errors. Their 
relation to the total yields of Kj^ which are established with considerably higher reliability up 
to SPS/Fermilab energies as well as their eventual comparison with strangeness production in 
nuclear collisions should therefore be critically reconsidered. 

3 The NA49 experiment, acceptance coverage and binning 

The basic features of the NA49 detectors have been described in detail in [[T]-[3l. The top 
view shown in Fig. |3] recalls the main components. 

The beam is a secondary hadron beam produced by 450 GeV/c primary protons imping- 
ing on a 10 cm long Be target. It is defined by a CEDAR Cerenkov counter, several scintillation 
counters (SI, S2, VO) and a set of high precision proportional chambers (BPDl-3). The hydro- 
gen target is placed in front of two superconducting Magnets (VTXl and VTX2). Four large 
volume Time Projection Chambers (VTPCl and VTPC2 inside the magnetic fields, MTPCL and 
MTPCR downstream of the magnets) provide for charged particle tracking and identification. 
A smaller Time Projection Chamber (GTPC) placed between the two magnets together with 
two Multiwire Proportional Chambers (VPCl and VPC2) in forward direction allows tracking 
in the high momentum region through the gaps between the principal track detectors. A Ring 
Calorimeter (RCal) closes the detector setup 18 m downstream of the target. 

The phase space region accessible to kaon detection is essentially only limited by the 
available number of 4.6 M inelastic events. It spans a range of transverse momenta between 0.05 
and 1.7 GeV/c for K+ and K and Feynman xp between and 0.5 for K~. For K"*" a limitation 
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Figure 3: NA49 detector layout and real tracks of a typical mean multiplicity p+p event. The 
open circles are the points registered in the TPC's, the dotted lines are the interpolated trajec- 
tories between the track segments and the extrapolations to the event vertex in the LH2 target. 
The beam and trigger definition counters are presented in the inset 



to Xi? < 0.4 is imposed by the constraints on particle identification discussed in Sect. |4]below. 

These kinematical regions are subdivided into bins in the Xplpx plane which vary ac- 
cording to the measured particle yields, effects of finite bin widths being corrected for in the 
evaluation of the inclusive cross sections (Sect. [5]). The resulting binning schemes are shown in 
Fig. m also indicating different ranges of the corresponding statistical errors. 
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Figure 4: Binning schemes mxp and pr for a) K"*" and b) K together with information on the 
statistical errors 
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4 Particle identification 



The identification of kaons by their ionization energy loss in the gas of the TPC detector 
system meets with specific problems if compared to pion |[T1 and baryon lO selection. This 
specificity has several reasons: 

- Corresponding to the momentum range of the NA49 data the ionization energy loss has 
to be determined in the region of the relativistic rise of the energy deposit, with the kaon 
energy loss positioned in between the one for baryons and for pions. 

- The relative distance in dE / dx between the different particle species is small and varies 
from only 4.5 to 7% for kaons with respect to protons and from 6.5 to 14% with respect 
to pions, over the xp range of the present data, with an rms width of the energy loss 
distributions of typically 3%. This creates an appreciable overlap problem over most of 
the phase space investigated. 

- High precision in the determination of the absolute position of the mean truncated en- 
ergy loss per particle species and of the corresponding widths is therefore mandatory. 

- The relative production yield of kaons is generally small as compared to pions, with 
K/tt ratios on the level of 5-30% for K+ and 5-20% for K~. In addition, for K+ the fast 
decrease of the K+/p ratio from typically 1 at x^r = to less than 5% atXF = 0.4 finally 
imposes a limit on the applicability of dE/dx identification towards high xp values. 
This general situation may be visualized by looking at a couple of typical dE/dx distri- 
butions for different xp regions as shown in Fig. [51 
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Figure 5: dE/dx distributions for K+ and K~ bins at xp = 0.05, pp = 0.4 GeV/c and xp = 0.25, 
Pt = 0.4 GeV/c superimposed with results of the fitted distributions 

As already described in [|2l| a considerable effort has been invested into the improved 
control of the analog response of the detector. Several aspects and results of this work, in partic- 
ular as far as kaon identification is concerned, will be discussed in the following sub-sections. 
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4.1 Non-Gaussian shape of the dE/dx distributions 

Due to the small K/vr and K/p ratios mentioned above, the precise description of the 
tails of the energy loss distributions of the dominant particle species becomes important. The 
extraction of kaon yields becomes indeed sensitive to small deviations in the upper tail of the 
proton and in the lower tail of the pion distributions for the extreme yield ratios mentioned 
above, as is also apparent from the examples shown in Fig.[5J Eventual asymmetries with respect 
to the generally assumed Gaussian shape of the energy loss distributions have therefore to be 
carefully investigated as they will influence both the fitted central position and the extracted 
yields of the kaons. A detailed study of the shape of the dE/dx distributions has therefore been 
performed both experimentally and by analytical calculation. 

By selecting long tracks in the NA49 TPC system which pass both through the VTPC 
and the MTPC detectors one may use the energy deposit in one of the TPC's to sharply select 
a specific particle type of high yield, for instance pions or protons. The dE/dx deposit in the 
other TPC will then allow a precise shape determination. An example is shown in Fig. |6]for the 
selection of pions atXF = 0.02 and px = 0.3 GeV/c in the VTPC. The corresponding distribution 
of the truncated mean for 90 samples in the MTPC is presented in Fig. [6^ together with a 
Gaussian fit. 




1.2 1.3 1.4 1.5 1.2 1.3 1.4 1.5 

dE/dx [MIP] dE/dx [MIP] 



Figure 6: a) Conventional Gaussian fit of the MTPC dE/dx distribution, for tracks with pion 
selection using the VTPC dE/dx; b) Ratio of data and fit function 



The small but very evident skewness of the truncated energy loss distribution is ex- 
pressed in Fig. by the ratio of the experimental data to a Gaussian fit. This ratio may be 
described by a cubic polynomial form with one normalization parameter Z, shown as the full 
line in Fig. [6b. 



(Data) /(Gaussian) ^ 1 + Z(/ - 3^), 



(3) 



where g is the distance from the mean of the dE/dx distribution, normalized to the rms of the 
Gaussian fit. 
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The parameter Z is related to the number of measured points, Ny, on each track, and the 
central dE/dx value by the relation 

with /3 and 7 experimentally determined to 0.5 [IJ and 0.4±0.2, respectively. Together with the 
relation: 



assuming = 7 which is a safe assumption regarding the sizeable error in the determination of 
7, 2:0 is obtained as 



zo = 0.215 ± 0.02 for the VTPC 

(7) 

^0 = 0.21 ±0.02 fortheMTPC. 

A Monte Carlo simulation based on the Photon Absorption Ionization (PAI) model ||3TI 
confirmed these results, demonstrating that the shape distortion is indeed a remnant of the basi- 
cally asymmetric Landau distribution of ionization energy loss. 

4.2 Position and width of the energy loss distributions 

Particle identification proceeds, in each defined bin of phase space, via a optimization 
procedure between the measured energy loss distributions and four single particle dE/dx dis- 
tributions of known shape but a priori unknown positions and widths for electrons, pions, kaons 
and protons, respectively. Due to the generally small fraction of electrons and their position in 
the density plateau of the energy loss function, and due to the known dependence of the dE/dx 
resolution on the dE/dx value for each particle species [[U, (Eq. (6]), the problem reduces in 
practice to the determination of eight quantities: three absolute positions of the energy loss of 
TT, K, p, one width parameter and four yield values which correspond to the particle cross sec- 
tions to be determined. If the fit of the pre-dominant particle species like pions and protons in 
general presents no problems, the situation is more critical for the kaons. Here it is in principle 
the central kaon position and the overall rms width of the dE/dx distributions which are liable 
to create systematic yield variations. In the ideal case, the detector response should reproduce 
exact scaling in the p/m variable as implied by the Bethe-Bloch function of ionization energy 
loss (BB), with p the lab momentum and m the particle mass. As shown in [[T]-[3l this scaling is 
fulfilled for pions and protons in the NA49 detector on the sub-percent level. The precision of 
the dE/dx fitting procedure allows for a quantification of the remnant deviations S with respect 
to the Bethe-Bloch parametrization as a function of xp and px 

6ixF,PT) = ^ixF,PT)- BB (8) 

in units of minimum ionization (MIP), where dE / dx is the mean truncated energy loss |[T1 . This 
is presented in Fig. |7] for the mean deviation of 7r+ and protons. 

The observed deviations are due to residual errors in the calibration of the detector re- 
sponse and in the transformation between the Bethe-Bloch parametrizations of the different 
gases used in the VTPC and MTPC detectors jSl. They stay in general below the level of 
±0.005. The fitted shifts of the kaon position, as characterized by their difference to the pion 



7 



J , , , I L 



0.5 1 1.5 

[GeV/c] 



Figure 7: Mean deviation {St,+ ^p), in units of minimum ionization, of 7r+ and proton dE/dx with 
respect to the Bethe-Bloch parametrization as a function of px for different values of a;i;' 



position 5k ~ ^ttj are shown in Fig. [8] as a function of xp and averaged over px, the error bars 
representing the rms deviation of the averages. 
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Figure 8: Mean deviations in units of minimum ionization of a) K+ and b) K with respect to 
the pion position {5f^± — (5^±) as a function of xp, averaged over px 



Evidently the measured positions fall well within the margin of ±0.005 in units of min- 
imum ionization as obtained for pions and protons. The similarity, within errors, between the 
results for and indicates systematic detector response effects as the principle source of 
the measured deviations. 

The fitted rms widths of the dE/dx distributions, characterized by their relative devia- 
tion from the calculated expectation value (Eq. [6] above), are shown in Fig.|9]as a function of 
Xp, after averaging overp^. 

The results show that the predicted widths are reproduced with an accuracy within a 
few percent of the expected values, with a slight systematic upwards trend as a function oi xp 
closely similar for K+ and K^. 
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Figure 9: Relative rms width function of xp for K"*" and K , averaged over px 



4.3 Estimation of systematic errors 

The dependence of the fitted kaon yields on the four parameters mentioned above, 
namely the positions of pions, kaons, protons, and the relative rms width of the fits, has been 
studied in detail. It appears that only two of these parameters are liable to produce noticeable 
systematic effects. These are the kaon position and the rms width. By enforcing a range of fixed 
values of these parameters, their influence on the extraction of kaon yields may be obtained. 
This is demonstrated in Fig. [10] for the dependence on kaon position and in Fig. [TT] for the de- 
pendence on the relative rms width, the error bars in each plot indicating the rms size of the pt 
dependence. 




Figure 10: Slope of the yield variation given in % per assumed kaon shift of 0.001 for K+ and 
as a function of x^, averaged overp-r 

Several aspects of this study are noteworthy: 

- As far as the influence of the kaon position uncertainty is concerned, and taking into 
account the size of the measured deviations from pions and protons and their rms fluc- 
tuation (see Fig. [8]) the related errors stay on the level of less than 1% up to = 0.2. 
Above this value the K+ yield reacts very critically on the fitted position. This is related 
to the proton yield which becomes rapidly overwhelming towards high xp. 

- Concerning the rms width the situation is somewhat more critical especially for K+. 
Here, allowing for a systematic error of about 0.5% in the fitted relative rms. Fig. [9l 
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Figure 11: Slope of the yield variation given in % per assumed change of a^^i of 0.005 for K"*" 
and K~ as a function of xp, averaged overp^ 



the corresponding yield error reaches values of about 2% at = 0.2 and about 10% at 
xp = 0.4. This is again measuring the influence of the large proton fraction. For on 
the other hand, the systematic error stays below the 2% level for the whole xp region 
investigated. 

The systematic errors estimated here have been included in the error estimation in Ta- 
ble E 



4.4 Fit stability and xp limit for kaon yield extraction 

The fitting procedure described above results in stable values for all eight parameters 
involved for xp values below about 0.25 for K+ and below 0.3 for K^. This is to be understood 
in the sense that the optimization procedure converges to a well-defined minimum in all 
variables with reasonable values for the ratio of over the degrees of freedom. For higher x p 
values the fits tend to become unstable in the sense that certain variables tend to "run away" 
into unphysical configurations. In the present case of extraction of kaon yields this concerns 
basically only the kaon position in the dE/dx variable and the rms width parameter of the 
energy loss distributions, as the pion and proton positions are always well constrained even in 
the critical regions of phase space. The problem is of course connected to the high sensitivity of 
the extracted kaon yield on these two parameters in relation to the small K/tt and K/p ratios as 
discussed in the preceding section. 

As the evolution of both the kaon position and the rms width with the phase space 
variables xp and shows no indication of any rapid variation up to the limits of fitting stability, 
and as indeed the geometrical configuration of the tracks in the TPC detectors shows a smooth 
and slow dependence on the track momenta in the regions concerned, it has been decided to 
extend the xp range up to 0.4 for K+ and to 0.5 for by imposing constraints on the two 
critical parameters. This is realized by constraining the kaon position to fixed values with respect 
to the pions, as indicated by the extrapolated lines in Fig.[8l and by also fixing the rms widths to 
the values following from the smooth extrapolation indicated in Fig. |9l The expected statistical 
error margins, allowing for reasonable values for the uncertainties in the quantities concerned, 
see Figs. [To] and [m have been added in quadrature to the statistical errors. 
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4.5 



Estimation of statistical errors 



It has been shown in f2\ that the estimation of the statistical error of the extracted particle 
yields has to take into account the dependence of the fit result on all parameters fitted via the 
covariance matrix. This means that the inverse square root of the predicted numbers of each 
particle species is only a first approximation to the relative statistical error. The fluctuations 
of the fitted particle positions discussed above and their contributions to the error of the yield 
parameters are intercorrelated with the particle ratios and with the relative distances of the 
energy deposits in the dE/dx variable. The method outlined in [2J has been applied to all 
extracted kaon yields and results in the statistical errors quoted in the data tables, Sect. [6]below. 
The ratio R^tat between the full statistical error and the inverse square root of the extracted 
yields is a sensitive indicator of the fluctuations inherent in the fitting method itself. It can vary 
drastically over phase space according to the correlation with the particle ratios and the relative 
positions with respect to the Bethe-Bloch function. This is visible in the distributions of the 
ratio i?stat defined above and shown in Fig. [121 for and in two different regions of x^. 




Figure 12: i?stat = osut/ for the bins xp < 0.2 (solid line) and 0.2< xp > 0.25 (dashed 
line); a) K+ and b) K" 



-Rstat is in general bigger for K+ than for due to the large p/p ratio. In both cases the 
forward bins in xp show a strong increase in i?stat which indicates the approach to the limit of 
stability of the fit procedure in particular for K"*". In the higher xp bins, xp = 0.3 and xp = 0.4 
the constraints imposed on some fit parameters. Sect. 14.41 limit of course also the range of the 
possible statistical fluctuations. Here, the problem has to be tracked by the evaluation of the 
corresponding systematic errors. 

5 Evaluation of invariant cross sections and corrections 

The experimental evaluation of the invariant cross section 

d^a 

f{xp,PT) = E{xp,PT)--j^{xp,PT) (9) 

follows the methods described in [1]. This includes the absolute normalization via the measured 
trigger cross section of 28.23 mb and the number of events originating from the liquid hydrogen 
target. The trigger is defined by a system of scintillation counters and proportional chambers on 
the incoming beam plus a downstream scintillator vetoing non-interacting beam particles. 



11 



5.1 Empty target correction 



Due to the small empty/full target ratio of 9% and the larger fraction of zero prong 
events in the empty target sample, the empty target contribution may be treated as a small 
correction as argued in [[TJ. This correction is, within the statistical errors, equal for and K~ 
and independent on and Xi;'. It is compatible with the one given for pions lUl and protons [|2l 
and is presented in Fig. [13] as a function oixp- 
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Figure 13: Empty target correction for K"*" and K as a function ofxp, averaged over 



5.2 Trigger bias correction 

This correction is necessitated by the interaction trigger which uses a small scintillator 
placed between the two magnets (S4 in Fig. [3]) in anti-coincidence with the beam signal. This 
trigger vetoes events with fast forward particles and thereby necessitates a trigger bias correction 
which can in principle depend both on particle type and on the kinematic variables. As described 
in detail in [jT|| the correction is quantified experimentally by increasing the diameter of the S4 
veto counter off-line and extrapolating the observed change in cross sections to diameter zero. 
For the case of kaons, the correction turns out to be within errors independent on px and similar 
for and K". Its xp dependence is shown in Fig. [141 




Xp Xp 

Figure 14: Trigger bias correction as a function ofxp for a) K+ and b) K~. The lines correspond 
to the parametrization of the correction 
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5.3 Re-interaction in the target 



This correction has been evaluated [IJ using the PYTHIA event generator. It is px inde- 
pendent within the available event statistics. The xp dependence is shown in Fig.fTSl 




0.2 0.4 



Figure 15: Target re-interaction correction as a function of xp 

5.4 Absorption in the detector material 

The correction for kaons interacting in the detector material downstream of the target is 
determined using the GEANT simulation of the NA49 detector, taking account of the K+ and 
inelastic cross sections in the mostly light nuclei (Air, Plastic foils. Ceramic rods). Due to 
the non-homogeneous material distribution the correction shows some structure both in pr and 
xp as presented in Fig.[T6l 




Xp 



Figure 16: Correction due to the absorption of produced kaons in the downstream detector 
material as a function of xp for two px values. The lines are shown to guide the eye 

5.5 Kaon weak decays 

Due to their decay length of about 30 m at the lowest lab momentum studied here, the 
weak decay of kaons necessitates corrections of up to 7% for for kaons produced in the lowest 
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measured xp range. Due to the high Q value of the decay channels and unlike the weak decay 
of pions, the decay products are not reconstructed to the primary vertex. This has been verified 
by detailed eyescans using identified kaons with visible decays inside the TPC system. The 
decay correction is therefore determined for those kaons which decay before having passed the 
necessary number of pad rows for reconstruction and identification. The resulting corrections 
are presented in Fig. [T71 
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Figure 17: Decay correction as a function of xp at different pr values. The lines are shown to 
guide the eye 



5.6 Binning correction 

The effect of finite bin sizes on the extracted inclusive cross sections is determined using 
the second derivatives of the xp or px distributions, as discussed in detail in [IJ. The associated 
corrections are within the statistical errors equal for K+ and K^. Examples are shown in Fig.fTSl 
as a function of at fixed pt and as a function of pr at fixed xp. 
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Figure 18: Binning correction, a) as a function of xp for px = 0.4 GeV/c and b) as a function 
of Pt for xp = 0.1. The crosses represent the corrections for fixed values of Sxp = 0.05 and 
SpT = 0. 1 GeV/c, respectively, and the open circles give the corrections for the used bin widths. 
The lines are shown to guide the eye 
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5.7 Systematic errors 



The systematic errors of the extracted cross sections are defined by the uncertainties of 
the normalization and correction procedures and by a contribution from particle identification as 
described in Sect. IH In particular the uncertainties due to the corrections may be well estimated 
from their distributions over all measured bins presented in Fig. \T9\ The corresponding error 
estimates are given in Table [B 
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Table 1 : Summary of systematic errors 



The linear sum of these estimations gives an upper limit of 4.5%, the quadratic sum 
an effective error of 2.2% for a;^ < 0.2. These values are close to the estimations obtained for 
pions [1] and baryons [2]. In the xp region above 0.25, however, the upper limit (quadratic sum) 
can reach 16.5%(12.2%) for K+ and 10.5%(6.4%) for K~. The spread of the corrections over 
all selected bins of phase space may be visualized in Fig. [19] which also gives the distribution 
of the sum of all corrections. 
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Figure 19: Distribution of correction for a) target re-interaction, b) trigger bias, c) absorption in 
detector material, d) kaon decay, e) empty target contribution, f) binning, g) total correction 
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Results on double differential cross sections 



6.1 Data tables 

The binning scheme described in Sect.[3]results in 158 data points each for K+ and K~. 
The corresponding cross sections are presented in Tables [2] and |3l 
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Table 2: Invariant cross section, f(xF,pT), in mb/(GeV^/c^) for K+ in p+p collisions at 
158 GeV/c beam momentum. The relative statistical errors, A/, are given in %, the transverse 
momentum px in GeV/c 
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Table 3: Invariant cross section, f{xp,pT), in mb/(GeV^/c'^) for K~ in p+p collisions at 
158 GeV/c beam momentum. The relative statistical errors. A/, are given in %, the transverse 
momentum pt in GeV/c 
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6.2 Interpolation scheme 



As in the preceding publications concerning p+p and p+C interactions |ni,(2l,l32l| a two- 
dimensional interpolation based on a multi-step recursive method using eyeball fits has been 
applied. The distribution of the differences of the measured points with respect to this interpo- 
lation, divided by the given statistical error should be Gaussian with mean zero and variance 
unity if the interpolation is bias-free and if the estimation of the statistical errors, see Sect. 14.51 
above, is correct. The corresponding distributions shown in Fig. [20] comply with this expecta- 
tion. 
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Figure 20: Difference A between the measured invariant cross sections and the corresponding 
interpolated values divided by the experimental uncertainty A/ for a) and b) 



As to first order the 8 first neighbours and to the second the 24 second neighbours to each 
data point contribute to the establishment of the interpolation, a reduction of the local statistical 
fluctuations of a factor of 3 to 4 may be expected. The authors therefore advise to use the data 
interpolation which is available under [33] for data comparison and analysis purposes. On the 
point-by-point level, the statistical error of the interpolated cross sections has been estimated as 
the mean value of the statistical errors of each measured point plus the 8 surrounding points in 
the PtIxf plane, divided by the (conservative) factor of 2. The systematic uncertainties are of 
course not touched by this procedure, in addition they are of mostly non-local origin. 

6.3 Dependence of invariant cross sections on xp and pt 

Shapes of the invariant cross sections as functions of pr and xp are shown in Figs. [2T] 
and [22] including the data interpolation presented above. In order to clearly demonstrate the 
shape evolution and to avoid overlap of plots and error bars, subsequent pt distributions have 
been multiplied by factors of 0.5 (Fig.[2T]). 

6.4 Rapidity and transverse mass distributions 

As in the preceding publications lfTl[2l[32l data and interpolation are also shown as func- 
tions of rapidity at fixed pt in Fig. [23] 



18 



CO 

CM 
> 

CD 
O 



E 



10-^ 



10 



■3 - 



10" 



10" 



: M 1 1 1 M 1 


1 1 M 1 M 1 : 


- a) 




r A. \ 


Xp = U.U 




(J.Ul 




0.025 - 




0.05 




0.075 - 




0.1 




\ 0.125 : 




A 0.15 


\ \ \ 
1 \ \ \ 

\ 

0.2 


\ \ \ \ \ \ 

\ \ \ \ \ V 
\ \ l\ \ \ " 

\ \ \' \ \ N 

\ \ \ \ \" 
\ \ \i \ 

\ \ \i \ - 

\ \ X \ " 


0.25 




- 0.3 


\ \ = 


0.4 

1 1 1 1 1 1 1 1 


1 \ \ : 

1 1 1 1 1 1 1 1 



0.5 



1 



1.5 




[GeV/c] 



[GeV/c] 



Figure 21: Double differential invariant cross section f{xF,PT) [mb/(GeV^/c^)] as a function 
of pt at fixed xp for a) K+ and b) K~ produced in p+p collisions at 158 GeV/c beam momen- 
tum. The distributions for different xp values are successively scaled down by 0.5 for better 
separation. The lines show the result of the data interpolation, Sect. 16.21 
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Figure 22: Double differential invariant cross section /{xp, Pt) [mb/(GeV^/c'^)] as a function of 
a;i? at fixed px for a) K+ and b) produced in p+p collisions at 158 GeV/c beam momentum. 
The lines show the result of the data interpolation, Sect. 16.21 
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Figure 23: Double differential invariant cross section f{xF,PT) [mb/(GeV^/c^] as a function of 
y at fixed px for a) K+ and b) produced in p+p collisions at 158 GeV/c beam momentum. 
The lines show the result of the data interpolation, Sect. 16.21 
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Figure 24: Invariant cross section as a function of — rriK for a) K+ and b) K produced at 
y = 0.0. The lines show the result of the data interpolation, Sect. 16.21 



Transverse mass distributions at xp = y = 0, with rriT = a/ m| + are presented in 
Fig. [24] again including the data interpolation. The corresponding dependence of the inverse 
slopes of these distributions on itlt — tuk is shown in Fig. [25] together with the results from the 
data interpolation. The local slope values are defined by three successive data points. 

K+ and K~ show a similar behaviour of the inverse slope parameters which fall from 
values around 180 MeV at low m^— m/< to a minimum of 150 MeV at m^— m/< ~ 0.05 GeV/c^ 
(Pt ~ 0.220 GeV/c), see also results from ISR (Fig. [74]). They then rise steadily towards higher 
TTiT and reach 180 MeV at pr ~ 0.6 (0.9) GeV/c and 200 MeV at ~ 1-0 (1.8) GeV/c for K+ 
and K , respectively. These trends resemble the ones observed for pions l[l|. 




Figure 25: Local inverse slope of the distribution as a function of — rriK for a) K+ and 
b) K~ . The lines correspond to the data interpolation, Sect. 16.21 



7 Particle ratios 

The present data on charged kaon production offer, together with the already available 
pion [[U and baryon ^ cross sections with similar phase space coverage and precision, a unique 
possibility to study particle ratios, in particular their evolution with transverse momentum and 
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X F ■ This section will therefore not only deal with K+/K but will also address K/tt and K/baryon 
ratios. 



7.1 K+/K ratios 



The ratio of the inclusive K+ and K cross sections, 



^K+K- = /(K+)//(K-) (10) 

is shown in Fig. |26] as a function of xp at fixed pt and in Fig. [27] as a function of px at fixed 
values of xp. 
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Figure 26: Rk+k- as a function of at fixed pt- The lines show the result of the data interpo- 
lation, Sect.O 



In each panel the corresponding ratio of the data interpolations. Sect. 16. 2[ is superim- 
posed to the data points as a solid line. The basic features of the data may be described as a 
steady increase of -Rk+k- over the available xp range by about a factor of three (Fig. [261) with 
some structure visible at certain xp and px values. The px dependence (Fig. [27]) reveals a more 
detailed evolution. The increase of -Rk+k- i^i the interval < pp < l.l GeV/c which amounts 
to about 60% at low xp flattens out in the xp range 0.1 - 0.2 to only 20% before it increases 
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again towards higher xp- This may be visualized in Fig. [28] where the ratios of the interpo- 
lated cross sections are shown as a function of pt for several xp values on an enlarged vertical 
scale. Fig. [28b gives an estimate of the statistical uncertainty of i?K+K- to be expected for data 
interpolation, characterized by the hatched area around the nominal values. 
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Figure 28: a) -Rr+k- for the data interpolation as a function of pr for different xp', b) Error 
bands expected for data interpolation 
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7.2 K/tt ratios 



The K/tt ratios shown here make use of the pion data and the corresponding interpolation 
published in [JJ . The ratios of the invariant inclusive cross sections 



i?K+.+ = /(K+)//(7r+) (11) 
i?K-.- = /(K-)//(7r-) (12) 

are presented in Figs. |29]to[34l 

Rfi+T,+ is shown in Fig. |29]as a function of px for fixed xp and in Fig. |30]as a function 
ofxp for fixed pt- 
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Figure 29: i?K+7r+ a function of pr at fixed xp- The lines show the result of the data interpo- 
lation, Sect. 16.21 




Figure 30: -RK+7r+ a function of xj? at fixed px- The lines show the result of the data interpo- 
lation, Sect. [O 



Here the salient features are the strong increase with pt which is rather independent 
on xp and reaches values of about 6 relative to low p-p at 1.7 GeV/c, and the rather small xp 
dependence with a slight increase at low pp and a comparable small decrease in the high pp 
region. These features are again shown in the composite Fig. [3T] where the pp dependence of 
i?K+7r+ from the interpolated data is plotted for the full range of values. 



25 



0.5 1 1.5 

[GeV/c] 

Figure 31: -RK+7r+ a function of px for different a; 

In Fig.|3T]the "cross-over" point at pt ~ 0.5-0.7 GeV/c where the full relative variation 
of i?K+7r+ with Xi? is on the level of only 20% of the ratio, and the practically parallel evolution 
of i?K+7r+ with Pt for different x f over a wide range of transverse momentum should be pointed 
out. 

R^-TT- is shown in Fig. [32] as a function of px at fixed x p and in Fig. |33]as a function of 
xp for fixed px- 
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Figure 32: -Rk-tt- as a function of pr at fixed xp- The lines show the result of the data interpo- 
lation, Sect. 16.21 

Also for -Rk-tt- a strong increase with px and the independence on xp for low px fol- 
lowed by a decrease with xp at high are evident. This is visualized in the composite Fig. [341 
where the px dependence for several xp values is plotted for the interpolated data values. 

Again a "cross-over" point in px with a practically complete xp independence, for 
i?K-7r- at Px ~ 0.3 GeV/c should be mentioned, together with the more pronounced decrease at 
higher Pt-. A general remark concerns the low px regions of Figs.[29l[3Tl[32]and [341 The rapid 
variation of the K/tt ratios below px ~ 0.2 GeV/c with some minima at px ~ 0.15 GeV/c are 
due to the structure of the vr^ and vr^ cross sections observed in this region [[U. This structure 
is more pronounced for 7r+ than for vr^ and has been explained by resonance decay |[Tl[34l. 
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Figure 33: -Rk-tt- a function of xj? at fixed px- The lines show the result of the data interpo- 
lation, Sect. [O 
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Figure 34: R^-^,- as a function of px for different xp 

7.3 K/baryon ratios 

The K/baryon ratios shown below use the new data on proton and anti-proton production 
published in The ratios of the invariant inclusive cross sections 

i?K+p = /(K+)//(p) (13) 
i?K-p = /(K-)//(p) (14) 

are presented in Figs. [35] to 1411 

i?K+p is shown in Fig. [35] as a function of pr for fixed xp and in Fig. [36|as a function of 
Xp for fixed px- 

Fig. [35] indicates a strong, rapidly decreasing K+ component at low pp and xp < 0.15, 
superimposed on an almost pt independent contribution which shows a marked decrease with 
increasing xp but also a slight increase with p-patxp > 0.2. This corresponds to the strong xp 
dependence at low pt in Fig. [36[ which flattens out rapidly with increasing px- The compos- 
ite Fig. [37] joins these trends using the ratio of the data interpolations as a function of pp for 
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Figure 35: -Rk+p a function of pt at fixed xp- The lines show the result of the data interpola- 
tion, Sect.O 
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Figure 36: i?K+p a function of Xi;' at fixed pt- The lines show the result of the data interpola- 
tion, Sect. [Q 

different xp values. 




Figure 37: a) -Rk+p for the data interpolation as a function of pt for different x^; b) Error bands 
expected for data interpolations 
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i?K+p seems to converge towards high pt to an a; ^ independent value of about 0.4 - 
0.5 as indicated in Fig. [37] by the dashed extrapolated lines for the different xp values. This 
is reminiscent of a similar behaviour for the p/(7r) ratio pointed out in [2J. As the point of 
convergence seems to lie close to p-r ~ 3 GeV/c it is tempting to use the available data at this 
transverse momentum from different a/s, although the detailed study of the s-dependence of 
i?K+p is outside the scope of this work. The analysis of the existing data at = 3 GeV/c and 
xf = from Serpukhov energy El via Fermilab [Sa to ISR ifTSUTTllIOHlIll. Fig. [M shows 
indeed consistency within errors with the value from the extrapolation shown above, indicating 
at the same time the very strong s-dependence of this particle ratio at high px- 
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Figure 38: s-dependence of -Rr+p at = 3 GeV/c and xj? = 0. The open circle corresponds to 
the NA49 extrapolation, Fig.ITT] 

It should be remarked here that the Fermilab data have been corrected for a systematic 
effect of 20% concerning the proton yields discussed in [2] and all ISR ratios by 10% to account 
for the expected amount of proton feed-down from strange baryons. 

i?K;-p is shown in Fig. [39] as a function of pt for fixed xf and in Fig. [40] as a function of 
Xp for fixed pt- 







Xp = 0.0 




Xp = 0.05 - 




Xp = 0.1 - 




Xp = 0.15 






















Xf = 0.2 




Xp = 0.25 - 




Xp = 0.3 - 










= 0.4 
















1 























0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 



p_^[GeV/c] P-r[GeV/c] P-|-[GeV/c] P-^IGeV/c] 

Figure 39: -Rk p ^s a function of px at fixed xp- The lines show the result of the data interpola- 
tion, Sect^^ 
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Figure 40: -Rr-p a function of at fixed pr- The lines show the result of the data interpola- 
tion, Sect.O 



As for -Rk+p the sizeable pr dependence at low xp flattens out at medium xp, 0.15 < 
xp < 0.25, and re-appears towards xp = 0.4. The x^? dependence, Fig.|40l is very different from 
the one of -Rk+p- There is no strong enhancement at low xp, R^-^ being rather independent on 
Xi;' up to Xi? ~ 0.3. Beyond this value the ratio increases rapidly towards values between 5 and 
6 at the maximum accessible xp. Fig. |4T] shows these trends using the ratio of the interpolated 
data as a function of px for different xp. 
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Figure 41 : a) -Rr-p for the data interpolation as a function of pr for different x^; b) Error bands 
expected for data interpolations 

The small structures described above are clearly visible, together with the strong in- 
crease at Xi? > 0.25 and a minimum dXpT values between 0.9 and 1.1 GeV/c. 
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Comparison to Fermilab data 



In a first step of data comparison, the NA49 data will be compared to the existing, 
double differential cross sections at neighbouring energies in order to control data consistency 
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with only small necessary corrections for s dependence. A wider range of comparisons ranging 
from kaon threshold to RHIC and collider energies will be performed in Sect. \W\ below. For 
the case of kaons all comparisons are facilitated by the absence of feed-down corrections from 
weak decays of strange particles. 



8.1 The Brenner et al. data fl3l 

This experiment which has shown a good agreement on the level of the double differ- 
ential cross sections for pions [IJ and baryons [2], offers 37 data points for K+ and 32 points 
for at the two beam momenta of 100 and 175 GeV/c. The average statistical errors of these 
data are unfortunately rather large for the kaon samples, with about 25% for K+ and 40% for 
K^. This is shown in the error distributions of Fig.|42l panels a) and d). Although the ^/s values 
at the two beam energies are, with 13.5 and 18.1 GeV, close to the NA49 energy, an upwards 
correction of 8% (12%) at the lower energy and a downwards correction of -2% (-5%) at the 
higher energy has been applied for and K", respectively, see Sect. [JO] for a more detailed 
discussion of s dependence. 




Figure 42: Statistical analysis of the difference between the measurements of |[T3l| and NA49 
for K+ (upper three panels) and K~ (lower three panels): a) and d) error of the difference of 
the measurements; b) and e) difference of the measurements; c) and f) difference divided by the 
error 



The statistical analysis of the differences between the Brenner et al. data and the interpo- 
lated NA49 results is presented in Fig.|42l Although the relative differences, dominated by the 
statistical errors of [[T3l . are very sizeable, see panels b) and e), the differences normalized to the 
given statistical errors, panels c) and f) show reasonable agreement between the two data sets. 
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in particular for K+ where the normalized differences are centered at A/cr = with the expected 
variance of unity. The show a positive offset of 0.3 standard deviations which corresponds 
to an average difference of 19%. 

A visualization of the Brenner data with respect to the interpolated NA49 results and 
their distribution in the xp and pr variables is given in Fig.l43l 




[GeV/c] Xp 

Figure 43: Comparison of invariant cross section between NA49 (lines) and measurements from 
[fT3l at 100 (full circles) and 175 GeV/c (open circles) for K+ as a function of a) px at fixed xp 
and b) Xi? at fixed pt, and for as a function of c) pt at fixed xp and d) Xi? at fixed pp. The 
data were successively divided by 4 for better separation 

Taking into account the comparison of all measured particle species for the two experi- 
ments OII2I it may be stated that a rather satisfactory overall agreement, within the limits of the 
respective systematic and statistical errors, has been demonstrated. 
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8.2 The Johnson et al. data ^ 



This experiment gives 40 data points for K+ and 50 points for within the range of the 
NA49 data obtained at 100, 200 and 400 GeV/c beam momentum. For comparison purposes the 
data have been corrected to 158 GeV/c beam momentum using the s-dependence established in 
Sect. [To] below. The distribution of the relative statistical errors is shown in Fig.l44lpanel a) for 
K"*" and panel d) for K^, with mean values of 12% and 9%, respectively. This is substantially 
below the errors of [JSil. 
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Figure 44: Statistical analysis of the difference between the measurements of [T7| and NA49 
for K+ (upper three panels) and (lower three panels): a) and d) error of the difference of 
the measurements; b) and e) difference of the measurements; c) and f) difference divided by the 
error 



The statistical analysis of the differences with respect to the interpolated cross sections 
of NA49 is also given in Fig.|44]in terms of the distribution of the relative difference A, panels 
b),e) and of the difference normalized to the statistical error A/ a, panels c) and f). Two main 
features are apparent from this comparison: an upwards shift of about 23% (10%) correspond- 
ing to 2.5 (1.0) standard deviations and large fluctuations corresponding to 1.7 (3.0) standard 
deviations for K+ and K^, respectively. As similar observations have been made for pions |[T1 
and baryons Q one may state that a general offset of 10 - 20% seems to be present which is 
compatible with the normalization uncertainty given in |[T4l . The fact that the proton data show 
a smaller offset might be connected with their x p coverage which is mostly at large negative x p 
(low lab momenta). On the other hand, the underestimation of the point by point fluctuations 
by a factor of 2 to 4 with respect to the claimed statistical errors, for all particle species, has to 
remain unresolved. 
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The phase space distribution of the data of lfT4l is shown in Fig.45 as a function of xp 
at fixed values of pr in comparison with the interpolated NA49 cross sections. 




Figure 45: Comparison of invariant cross section between NA49 (lines) and measurements from 
[fT4]| at 100 (full circles), 200 (open circles) and 400 GeV/c (full triangles) as a function of xp 
at fixed pt for a) and b) K^. The data were successively divided by 3 for better separation 



8.3 The Antreasyan et al. data lHHl 

It is only the low-p^ part of this experiment which can be compared to the NA49 data, 
at xp close to 0. Due to the fact that the spectrometer of ffT2ll was set to a constant lab angle 
for all beam energies and particle species, the given cross sections have to be compared at their 
proper Xi? values as given in Table H see also the corresponding arguments in tlj. 



PT [GeV/c] 


Pbeam [GeV/c] 

[GeV] 


200 
19.3 


300 
23.7 


400 
27.3 




Xf 


-0.0054 


-0.011 


-0.020 


0.77 


Rk+ 


0.826±0.12 


1.026±0.16 


1.110±0.20 




Rk- 


0.966±0.12 


1.217±0.18 


1.164±0.18 


1.54 


xf 
Rk+ 


0.0302 
0.796±0.05 
0.791±0.06 


-0.031 
1.080±0.08 
1.240±0.06 


-0.020 
1.260±0.12 
1.616±0.14 



Table 4: Offset in at different y/s and pp- The cross section ratio Rf^± between the data 
from [12] andNA49. 

The cross section ratios Rj^+ and R^^- are shown in Fig. |46]as a function of y/s at fixed 
Pt, together with the s-dependence extracted in Sect.[10]below from data at a; i? = at Serpukhov 
energy mUl and ISR energy [12111221 . 

Evidently the data [fT2l comply, within their sizeable statistical errors, with the s- 
dependence as established by the other experiments. However, three of the four points at 
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Figure 46: The cross section ratios between the data from [fT2l and NA49 as a function of y/s 
for two values of px for a) and b) K . In both of the panels the NA49 point is indicated 
with full triangle. The full and dashed lines represent the result of the s-dependence at = 
established in Sect.[10]below atpT = 0.77 and 1.54 GeV/c, respectively 



200 GeV/c beam momentum are low by about two standard deviations. This would, by us- 
ing the data [12J alone to establish the s-dependence, lead to a large underestimation of the 
kaon yields at lower s. See also the discussion in for baryons. 

8.4 Comparison of particle ratios 

As systematic effects tend in general to be reduced in particle ratios, it is interesting to 
also look at the consistency of the corresponding ratios from lfT2] - [T4l with the NA49 data shown 
in Sect. |7] of this paper. This is shown in Fig.|47]for -Rr+k- ratios, in Fig.|48]for i?K±7r± 



40 



20 



[12] 
* 200 GeV/c 
A 300 GeV/c 
» 400 GeV/c 

0.4 



0.3 



0.2 



0.0 



_l I I I I I I L_ 



60 



40 



20 



- b) 


1 1 1 1 1 1 1 1 1 1 

[13] 
• 100 GeV/c 
o 1 75 GeV/c 


1 1 1 1 1 1 1 1 1 

[14] 
■ 100 GeV/c 
° 200 GeV/c 


- Pt- 


1 .0 GeV/c 




0.75 




£ 








0.3 a * ? ' 


n qe; 





, * ; 


1 1 1 1 1 rr-rri , , , , i , , , , i , , 


) 


0.1 0.2 


0.3 0.4 



0.5 1 1.5 

[GeV/c] 



Figure 47: Comparison of Rk+k- between [fT2l (triangles), (13] (circles), [fT4l (squares) and 
NA49 (lines) as a function of a) and b) pr- The data were successively shifted by 10 for 
better separation 
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Figure 48: Comparison between (T2\ (triangles), (13] (circles), [fT4l (squares) and NA49 (lines) 
of -RK+7r+ ^ function of a) and h) xp and -Rk tt- ^ function of c) and d) xp- The data 
were successively shifted by 0.2 for better separation 



Fig.|49]for K/baryon ratios. 



8.5 Conclusion from data comparison at Fermilab energies 

In conclusion of the detailed comparisons in the Fermilab/SPS energy range shown 
above it may be stated that a mutually consistent picture for kaon production from several 
independent experiments has been established, with the exception of some offsets in the ab- 
solute cross section especially for [fT2ll and lfT4l . These offsets tend to cancel in the particle 
ratios -Rr+k- for both lfT2]| and lfT4]| . The ratios Rkw and -Rkp are consistent for [131 and lfT4l 
within their statistical uncertainties, whereas for lfT2ll the systematic effects discussed in ^ for 
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Figure 49: Comparison between [T2^ (triangles), (13] (circles), [fT4l (squares) and NA49 (lines) 
of -Rk+p a function of a) px and b) a;^? and -Rk^p ^ function of c) px and d) x^?. The data 
were successively shifted by 1.2 for -Rk+p ^i^d by 6 for -Rk-p for better separation 



baryons and in Sect. l8.3l for kaons persist for i?K7r and _Rkp- What is also important to note is the 
apparent absence of systematic deviations as a function of kinematic variables xp and px- This 
lends, as none of the existing experiments has on its own sufficient phase space coverage, some 
confidence to the establishment of pr integrated and total yields from the NA49 measurements 
alone, as discussed below. 



9 Integrated data 

In a first step the data interpolation. Sect. 16. 2[ will be used to perform an integration 
over transverse momentum. In a second step the total charged kaon yields will be determined. 
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These can be used, in conjunction with the total pion and baryon yields published before [[T1I3 
to control the total charged multiplicity with respect to the precision data from bubble chamber 
experiments. 

9.1 Pt integrated distributions 

The Pt integrated non-invariant and invariant kaon yields are defined by: 

dXF CTinel 2 J E 

F = j fdvl (15) 
dn TT /■ , , 2 

dy (Tinel J 

with f = E ■ d^a/dp^, the invariant double differential cross section. The integrations are per- 
formed numerically using the two-dimensional data interpolation (Sect. 16.21) which is available 
in steps of 0.05 GeV/c in transverse momentum. 
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Table 5: pt integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum {pt) [GeV/c], mean transverse momentum squared (p^) [(GeV/c)^] as 
a function of x^, as well as density distribution dn/dy as a function of y for K+ and K^. The 
statistical uncertainty A for each quantity is given in % as an upper limit considering the full 
statistical error of each measured PtIxf bin 

The statistical uncertainties of the integrated quantities given in Table [5] are upper lim- 
its obtained by using the full statistical fluctuations over the measured bins. As such they are 
equivalent, for the kaon yields, to the statistical error of the total number of kaons contained in 
each Xp bin. 

The resulting distributions are shown in Fig. [50] for K+ and as a function of and 
y. The relative statistical errors of all quantities are generally below the percent level. They 
increase towards the high end of the available xp region essentially defined by the available 
event number and, especially for K+, by limits concerning particle identification (Sect.H]). The 
K+/K~ ratio, (pr) and {p\) for kaons as a function of are presented in Fig.[5Tk-c. Fig. [5Tb 
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shows the mean transverse momentum of all measured particle species in a single panel in 
order to allow a general overview of the interesting evolution of this quantity with xp which 
demonstrates that (p^) is equal to within 0.05 GeV/c for all particles at xp ~ 0.3 - 0.4. 




y 



Figure 50: Integrated distributions of K+ and K produced in p+p interactions at 158 GeV/c: 
a) density distribution dn/dxp asa function of xp', b) invariant cross section F as a function of 
xp\c) density distribution dn/dy as a function of y 



+ 




0.2 0.4 0.6 0.8 



Figure 51: a) K+/K^ ratio, b) mean pp, and c) mean p'^ as a function of xp for K+ and 
produced in p+p interactions at 158 GeV/c; d) mean pp for 7r+, tt^, K+, K^, p, p on an enlarged 
vertical scale 



9.2 Comparison to other data 

As in Sect. [HI a first stage of the comparison is limited to the SPS/Fermilab energy range 
where only two experiments provide integrated cross sections. The data of Brenner et al. |fT3l 
are obtained from a limited set of double differential cross sections, using basically exponential 
fits to the measured points. The resulting invariant cross sections F(xp) are shown in Fig.[52lin 
comparison to the NA49 data. 

As already remarked for the case of pions and protons, very sizeable deviations are 
visible in the distributions of Fig. [52k. which are quantified by the ratio of the two measurements 
shown in Fig. [52b . If the relative differences in F were limited to about ±40% for pions and 
protons [[BO, the factors are even bigger for kaons, with a mean deviation of about 50%. This 
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Figure 52: a) Comparison of pt integrated invariant cross section F as a function ofxp for K"*" 
and measured by [13] to NA49 results (represented as lines). The data for are multiplied 
by 0.1; b) Ratio R between measurements of [|T3l and NA49 results. The mean ratios for K+ 
and are presented with dashed lines 

again demonstrates the danger of using oversimplified algebraic parametrizations of double 
differential data which comply with the NA49 measurements on the point-by-point level within 
their statistical errors (Sect. 18. ll) . 

The EHS experiment [35] at the CERN SPS, using a 400 GeV/c proton beam, offers 
integrated data which are directly comparable in all quantities defined in Eq. [151 In view of the 
s-dependence which is enhanced at low in the quantity dn/dxp [11]I2J and of the important 
shape change to be expected in the rapidity distributions, only the invariant integrated cross 
section F is plotted in Fig.[53]in comparison to the NA49 data. 




Figure 53: Comparison of pr integrated invariant cross section F as a function of for a) K+ 
and b) measured by [i35il to NA49 results (represented as dashed lines) 

Some remarks are in place here. The EHS K+ data show an enhancement at low xp oi 
about 35% which is substantially above the expected s-dependence, see also the discussion of 
the kaon data in Sect. [131 After a local deviation from a smooth x p dependence at ~ 0. 1 5 the 
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distribution cuts, however, below the NA49 data in the region 0.175 < xp < 0.45. This decrease 
cannot be explained by any known s-dependence. For K the situation is qualitatively similar. 
Here at s^r = an enhancement of 41% is observed, with an xp dependence which smoothly 
approaches the NA49 data to become equal to these cross sections within errors at xp > 0.22. 
Again such behaviour contradicts the expected s-dependence. A possible explanation might be 
contained in the mean p"^ data shown in Fig. [54l If the results on (pf.) agree at xp = within 
the respective errors, the EHS data deviate rapidly upwards from the NA49 measurements with 
increasing xp. For K+ the instability in the cross sections at Xi? = 0.15 is seen as a break in the 
Xp dependence of at the same xp value. The xp dependence then flattens in the region 
0.2 < Xp < 0.45 which corresponds to the depletion of the cross section, rising again steeply to 
very large values at xp beyond the range accessible to NA49. A similar behaviour is observed 
for where (p^) shows reasonable consistency up to Xi;' ~ 0.2 with a slight increase over 
the NA49 data which is however inconsistent with the s-dependence in Sect. 110.6.41 Above 
Xp 0.2, however, there is again a strong almost linear increase of with xp with values in 
excess of 0.8 (GeV/c)^ in the high xp region. One may speculate that both the behaviour of the 
invariant cross sections and the one of (p^) are of the same origin if one assumes that there are 
detection losses for kaons with increasing xp and at transverse momenta below the mean value. 
This would reduce the observed cross sections and enhance the mean p"^. 




Figure 54: Comparison {p"^) as a function of for a) K+ and b) measured by [[351 to NA49 
results (represented as dashed lines) 

In conclusion of the comparisons with the EHS experiment which have been carried out 
with some precision for pions [IJ, baryons [2J and here for kaons, a somewhat unsatisfactory 
and partially inconsistent picture emerges. In general it may be stated that sizeable relative 
differences, even after taking into account possible s-dependences, emerge at a level of typically 
±10 - 30% which cannot be explained by a common factor like normalization uncertainties. In 
addition there seems to be a general tendency of unphysical behaviour in the EHS data for x p 
values above about 0.2 both in the cross sections and, more extremely, for the behaviour of 
mean p"^. 

9.3 Total kaon yields and mean charged multiplicity 

For the x p integration of the dn/dxp distributions presented in Table |5] an exponential 
extrapolation into the unmeasured region at high xj? has been used. This is well justified by the 
shape of the distributions within the measured region and by the fact that only 4% (0.3%) of the 
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total yields are beyond the experimental limits for K+ and K , respectively. The resulting total 
kaon yields are: 

(nK+) = 0.2267 
(riK-) = 0.1303 

= 1-740 (16) 

2 

The statistical errors of these yields may be estimated by the total number of kaons 
extracted from the 4.8M events of this experiment. These are 260k for K+ and 170k for K^. 
From these numbers follows, including the additional statistical errors from particle identifica- 
tion, Sect. l4.5I Fig.[T2l an error of 0.27% for K+ and 0.28% for K~ which is about one order of 
magnitude below the smallest estimated systematic error (Table [!])• 

These numbers, together with the results for pions [IJ and baryons {T\ can be used 
to establish the mean charged multiplicity as it results from this experiment. The respective 
numbers are given in Table [6] below. 





positives 


negatives 


total 




3.018 


2.360 


5.378 




0.227 


0.130 


0.357 


(^p) 


1.162 


0.039 


1.201 


in) 


4.407 


2.529 


6.936 



Table 6: Mean multiplicities of charged particles 

In order to establish the total charged multiplicity and to be able to compare to the results 
from bubble chamber work where the charged hyperons are included as on- vertex tracks, an 
estimation of S+ and yields has to be performed. Several measurements of S+, and 
S° are available in the energy range 3 < y/s < 27 GeV, all with rather big relative statistical 
errors of typically 15 to 50%. For the present purpose where the charged hyperons constitute 
a correction of about 1%, this is nevertheless acceptable since all results stem from bubble 
chamber experiments with small systematic uncertainties. Three quantities are interesting and 
necessary for the present comparison: 

1 . the S°/A ratio 

2. the ratio 

3. the ratio (S" -i- S° -i- S+)/A 

The S°/A ratio has been obtained by 5 experiments If36] - l40]| with values between 0.1 
and 0.74 with an average of 0.4. This value may be used to obtain the ratio (S^ -i- S° -i- S+)/A 
[|4Tti43]l which varies between 0.83 and 1.09 with an average of 0.99. The ratio S+ZS" Il39ll411- 
|43l|45l| shows a variation from 2 to 5.2 with an average of 3.3. 

Adopting the average values for the ratios (1) and (3) the combined yield S"*" -i- E~ may 
be obtained at = 17.2 GeV by interpolating the well-established total yield of A [I39l - l44ll46l - 
|59] to {nj\) = 0.12 per inelastic event at this energy. This results in a contribution of 0.07 per 
inelastic event from charged hyperons and gives a total charged multiplicity 

(rich) = 7.01 (17) 
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from this experiment. This multiplicity may be compared to the existing measurements essen- 
tially from Bubble Chamber experiments taken from |60] and presented in Fig. |551 




Figure 55: {rich) as a function of beam momentum pbeam- The NA49 measurement is indicated 
with an open circle 

The full line in Fig.[55]represents a hand interpolation of the measurements in the range 
from 50 to 300 GeV/c beam momentum. It coincides incidentally, at ^/s = 17.2 GeV, with the 
parametrization 

(rieh) = -4.8 + 10/v^ + 2.01ns (18) 

given by [6T| which predicts 

(nch) = 7.15 (19) 

The relative deviation of the summed integrated yields given above from this value cor- 
responds to -2%. It is certainly governed by the systematic uncertainties of the dominant pion 
and proton yields for which the systematic error estimation [1, 2] gave 4.8% (5%) for the lin- 
ear sum and 2% (2.5%) for the more optimistic quadratic sum of the contributions,respectively. 
Allowing for a typical error of about 1 % of the bubble chamber data, it may be stated that the 
observed deviation is within the error estimate for the NA49 data. 

At this point it is indicated to also check the charge balance of the NA49 results where 
the difference between positive and negative particle yields should give two units from charge 
conservation. Using the total charged hyperon yield estimated above and the average 
ratio of 3.3 the following yields are obtained: 



(ns+) = 0.054 
(ns-) = 0.016 
(ripo,) = 4.461 
(nneg) = 2.545 

(^pos) - (?^neg) = 1-916 



(20) 
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This means that the charge balance is off by 0.08 units or about 4% of its nominal value. 
In order to put this number into perspective it should be realized that a systematic downwards 
deviation of the 7r+ yield by 1.5% accompanied by an upward shift of the tt^ yield by the same 
relative amount is sufficient to explain this imbalance. Therefore it may be stated that also the 
charge conservation of the NA49 results is established within the stated systematic errors. 

10 A new evaluation of s dependence 

The new set of kaon data presented and discussed above has been used, in connection 
with existing data at other cms energies, to re-assess the experimental situation as far as the 
s-dependence, in particular also for integrated yields, is concerned. It is indeed rather surprising 
that the very first attempt in this direction by Rossi et al. [4] which dates from 1975, is still 
being used as a reference for rather far-reaching conclusions with respect to kaon production in 
heavy ion interactions [162]| . This is especially concerning the admitted systematic uncertainties 
which are given in dH as only 15% for their estimated total yields. In view of the rather sparse 
phase space coverage of most of the preceding data sets, see Sect.|2]and Fig.[T]above, it is in fact 
for most cms energies quite difficult to establish integrated yields with defendable reliability. 
In this context it is interesting to also look at the available data on K^, production which, com- 
ing for the s-range up to medium ISR energies exclusively from bubble chamber experiments, 
have well defined systematic errors in particular for integrated yields, notwithstanding their in 
general rather limited statistical significance. Here, the relation between charged and neutral 
kaon production deserves special attention as it is directly sensitive to the respective produc- 
tion mechanisms. In the following section, five energy ranges from Cosmotron up to RHIC and 
collider energies will be inspected in an attempt at establishing some coherence with respect to 
s-dependence. 

10.1 The K+ data of Hogan et al. [SJ and Reed et al. [6| at = 2.9 GeV 

These early experiments at the Princeton-Penn (PPA) and BNL Cosmotron accelerators 
use a range of beam momenta from 3.2 to 3.9 GeV/c with a common point at about 3.7 GeV/c. 
The data at this energy have been used in order to establish a maximum of combined phase 
space coverage in the ranges < xp < 0.4 and < px < 0.6 GeV/c, see Fig. [T^. It should 
be mentioned here that the definition of xf (Eq. |2l) has been used throughout although it pro- 
gressively limits the available xp range at low interaction energies due to energy-momentum 
conservation, see [8] for a detailed discussion. At y/s = 3 GeV this means rather sharp cut-offs 
in production cross section towards x^? ~ 0.5 and px ~ 0.7 GeV/c. Within these limits, reason- 
able inter- or extra-extrapolation may be performed in order to establish approximate px and xp 
dependences. It should be stressed that throughout this paper no arithmetic parametrizations of 
Xp or pt distributions have been used as those would introduce large systematic biases which 
are difficult to control. Instead, two-dimensional interpolation by multi-step eyeball fits, as dis- 
cussed in Sect. 16.21 above, have been applied. Two examples of this procedure are shown in 
Fig. [56] for selected xp and pt values, where the available, interpolated or slightly extrapolated 
data points are indicated. The resulting interpolation of cross sections over the complete Xi;' and 
Pt ranges is presented in Fig.[57l 

The interpolation shown in Fig. [57] may be pp integrated in order to obtain the F, 
dn/dxp and (pp) dependences shown in Fig. [58] 

In a second step the integration over xp may be performed resulting in an average 
K+ multiplicity of (nK+) = 0.00481. This value is 8.3% (6.0%) higher than the multiplicities 
(nK+) = 0.00441±17% and {n^+) = 0.00452±23% obtained by §^ and respectively. These 
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Figure 56: / as a function of px for a.) xp = and h) xp = 0.2. The data points from Hogan et 
al. O and Reed et al. [6J are given together with the data interpolation (full lines) 



Figure 57: Interpolated invariant cross sections as a function of a^i? for fixed values of px 

groups imposed isotropy (S-wave decay) in the cms system in order to be able to carry out the 
data integration. There is also a bubble chamber experiment from the BNL Cosmotron at the 
same beam momentum ll63l which gives (nK+) = 0.00462±19% for the multiplicity which 
is only 4% lower than the result obtained above. In conclusion a statistically consistent K+ yield 
from 3 independent experiments may be claimed at y/s = 3 GeV/c which is about 55% above 
the one elaborated in [4J . 

The bubble chamber experiment [l63ll also gives the K° multiplicity as (n^o) = 0.00165. 
With the usual assumption {n^o ) = 0.5 (nKo) this corresponds to {n^o ) = 0.000824. The ratio 
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Figure 58: integrated a) F, b) dn/dxp and c) {px) distributions as a function of xp- The 
results obtained at ^ = 17.2 GeV (dashed lines) are also shown for comparison 



is therefore, with {nj^+) from [163]| as given above, 2.8 which is substantially above the value 
-Rk|k± = 1 expected from isospin invariance. Inspecting the K° and data of [3T] and [|64l . 
_RfjOfj± is determined to 1.4 at ^/s = 3.5 GeV and 1.27 at ^/s = 4 GeV. This indicates a steep 
deviation from isospin invariance in kaon production as the threshold is approached from above, 
as shown in Fig. 
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Figure 59: Ratio -RkOk± between the average charged kaon and yields as a function of ^/s. 
The threshold of kaon production is indicated at about ^/s ~ 2.5 GeV 



Evidently -Rk^k± approaches unity rather quickly with increasing energy so that 
-Rro K± = 1 iTiay be assumed within a few percent error margin at y/s > 5 GeV, see Sect. [TT] 
below for a more detailed discussion. 

It is also interesting to compare the differential data of f5] and directly to the NA49 
data. The ratio of the invariant inclusive cross sections. 



f{xF,PT, ^/s = 3 GeV) 



(22) 



is shown in Fig. 

of pt ■ 



f{xF,PT,V^=n.2 GeV) 
as a function of pr at constant xp and as a function of Xi? at constant values 



46 



« 0.2| 




Figure 60: Ratio Rs as a function of a) pr at fixed xp and h) xp at fixed 

Evidently the total yield ratio of 0.021 does not translate into a common suppression 
factor for the differential distributions but the local cross section ratios show a strong and com- 
plex dependence on the kinematical variables. If the complete suppression of K+ production for 
Pt ^ 0.7 GeV/c and xp ^ 0.5 is a trivial consequence of energy-momentum conservation, the 
local structures as for instance the maximum at ~ 0.3 and low pr are a consequence of the 
evolution of different production mechanisms with increasing interaction energy. 



10.2 Data in the PS/AGS energy range 

In this subsection data in a range from 12.5 to 24 GeV/c beam momentum are grouped 
together, again in an effort to consolidate the available information and to quantify the consis- 
tency of the different data sets. This concerns the double differential data by Akerlof et al. [9] at 
12.5 GeV/c beam momentum, of Dekkers et al. ifTOl at 18.8 and 23.1 GeV/c, and the extensive 
data sets of the CERN/Rome group, AUaby et al. at 14.2, 19.2 and 24 GeV/c beam mo- 
mentum. The data sets from all these groups have been tabulated conveniently by Diddens and 
Schliipmann in Landoldt- Bomstein [651. As the overview of Fig.dt shows, there is a fair cov- 
erage of phase space and some mutual overlap, unfortunately again [1,2] with the exception of 
the low xp region, xp < 0.1-0.15, at all pt- In a first step, the AUaby et al. data [|7l[8|| are trans- 
formed to the standard xp values following Eq. [2|and interpolated using the two-dimensional, 
multistep eyeball method described in sect. 16.21 An extrapolation into the non-measured phase 
space areas is then attempted in order to allow the establishment of integrated yields. The situ- 
ation may be judged from Fig. [611 where the interpolated/extrapolated cross sections are shown 
as a function of xp at fixed values of Pt for K+, Fig. [6Tb. and K^, Fig. [6Tb. Here the regions 
with available measurements from [[7l[8l are marked by the hatched areas. 

Clearly, the above remark concerning the problems with data extrapolation is well in 
place here, especially for the higher pt regions. However, at least towards low pr there is not 
much freedom of choice, as well as for the px region < < 1 GeV/c towards xp = 0. Up 
to ~ 1 GeV/c it is hard to imagine an extrapolation which would be off by more than, say, 
10-20% from the lines shown at Xi? = 0. It is also clear that the increasing error margin towards 
higher px will not contribute too much to the integrated cross sections. The ratio 
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^ _ f{xF,pT,V~s = 6.SGoY) 
' /(xi.,pT,yi= 17.2GeV) 

is shown in Fig. |62]for and K+ as a function of px and xp- 
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Figure 62: Rg a) as a function of pt at fixed xp and b) as a function of for fixed px for K 
and c) as a function of px at fixed x p and d) as a function of x^? for fixed px for K+ 
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As already visible for the low energy data, Fig. [60] above, very clear structures are appar- 
ent with local maxima at pt < 0.5 GeV/c and xp around 0.2 - 0.4. Evidently the s-dependence 
varies over phase space by more than an order of magnitude. The apparent s-independence, 
within 5%, of the K+ cross sections in the region oixp around 0.3 for p-p from up to 0.6 GeV/c, 
should however be seen as an indicator of systematic problems. The fact that the data of f7ll8| 
are apparently over-estimating the K+ yields in this region can be shown in comparison with 
data from other experiments [|9llT0l. Although these data do not offer enough coverage to permit 
a complete interpolation, they may be used to bring out local mutual inconsistencies between 
the experimental results. As the data of Dekkers et al. [10] have been obtained at 18.8 and 
23.1 GeV/c beam momenta, at only two fixed lab angles of and 100 mrad, also the AUaby 
et al. data at 19.2 GeV/c beam momentum, which offer much inferior phase space coverage, 
have been interpolated in order to permit direct comparison for a maximum of data points. The 
results are presented in Figs. [63] and [64] 
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Figure 63: comparison Dekkers [lOJ and AUaby a) K+, px = GeV/c as a function of 
xp- Full lines [[71[8l extrapolation at 19.2 and 24 GeV/c beam momentum, data points from 
[[T0]| at 18.8 and 23.1 GeV/c b) Dekkers data and AUaby interpolation at 18.8 (19.2) GeV/c 
and 100 mrad lab angle as a function of xp- Full lines AUaby interpolation, data points from 
Dekkers. The data at 23. 1 GeV/c and lines at 24 GeV/c are multiplied by 0. 1 for better separation 

If the degree data (Fig. [63k ) show all the Dekkers points below the AUaby extrapola- 
tion, with a mean relative difference of 20%, the values at higher pr. Fig. [63b . are all far below 
the interpolation by a factor of about 2. The same comparison for is shown in Fig. [64] 

Also for the Dekkers data dXpx = GeV/c are below the AUaby data, here by 30%. 
This might indicate a general offset between the two data sets of 20-30% which does not seem to 
be excluded by the systematic uncertainties given for the respective experiments. In contrast to 
the situation for K"*" the data interpolation at 19.2 GeV/c at higher px is bracketing the Dekkers 
data for K~ such that the mean deviation over the given px scale tends to be small. Fig. l64b . 

A further possibility of controlling s-dependence is given by the data of Akerlof et al. |9l 
which were obtained with a 12.5 GeV/c proton beam at the Argonne ZGS. Although only 7 
points for and 17 points for have been measured, the s-dependence between these data 
is revealing if compared to the 19.2 and 24 GeV/c data of AUaby et al. Starting with K^, the 
Akerlof data allow comparison at fixed pt of 0.632 GeV/c in the xp range from 0.12 to 0.32, 
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Xp = 0.259 
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Figure 64: comparison Dekkers [lOJ and Allaby a) K^, pr = GeV/c as a function of 
Xp- Full lines [|71[8l extrapolation at 19.2 and 24 GeV/c beam momentum, data points from 
[fTOl at 18.8 and 23.1 GeV/c b) Dekkers data and Allaby interpolation at 18.8 (19.2) GeV/c 
and 100 mrad lab angle as a function of xp. Full lines Allaby interpolation, data points from 
Dekkers. The data at 23. 1 GeV/c and lines at 24 GeV/c are multiplied by 0. 1 for better separation 



and at fixed xp = 0-24 for pt between 0.55 and 1.14 GeV/c. Fig. [65] shows the s -dependence 
for fixed xp (panel a) and fixed pp (panel b) including the data from Allaby and NA49. 



O 

CM 
> 

O 

E 10-^ 



10"' 



: 0.632 GeV/c 



K" 




0.2 



Xp = 0.32 



o[9] 

i [8], 19.2GeV/c 
A [8], 24 GeV/c . 
• NA49 



10 20 

\Js [GeV] 



10" 



10' 



10' 



: b) 




. Xp = 0.24 


P^IGeV/c] . 




0.55 _ 




0.8 







5 10 


20 



\ls [GeV] 



Figure 65: K comparison as a function of y/s, a) pp = 0.632 GeV/c for Xi? = 0.12, 0.2 and 0.32 
b) Xp = 0.24 for pt = 0.55, 0.8 and 1.0 GeV/c 



For all xp/pp combinations, a smooth s-dependence between the four data sets is 
observed. A different picture emerges for K+, Fig. [66] where only 3 points in y/s at 
Pp = 0.632 GeV/c and xp = 0.24 are available. 

As already apparent from Fig. [62] the Allaby et al data at pp = 0.632 GeV/c, = 0.2 
(Fig. [66b ) andxp = 0.24, pp = 0.55 GeV/c (Fig. [6^5) are on the same level as the NA49 data 
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Figure 66: comparison as a function of y/s a) px = 0.632 GeV/c for xp = 0.2 h) xp = 0.24 
forpT = 0.55, 1.18 and 1.26 GeV/c 



for 24 GeV/c beam momentum, and are even higher for 19.2 GeV/c. Compared to this the 
Akerlof data show the expected decrease at y/s = 5 GeV. The resulting s-dependence looks 
definitely unphysical indicating an excess of the order of 60% in the K+ yields of AUaby et al. 
A similar problem is present in the AUaby et al. data at 14.25 GeV/c beam momentum which 
only exist for K+ at a lab angle of 12 mrad, thus covering the low px region from 0.04 to 
0. 1 1 GeV/c for 0.25 < xp < 0.6. Those data may be compared to the interpolation at 24 GeV/c 
beam momentum. The cross section ratio /(14.25 GeV/c)//(24 GeV/c) is shown in Fig. [67] as 
a function of x^. 
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Figure 67: /(14.25 GeV/c)//(24 GeV/c) as a function of x^. 



Evidently there is a very small s-dependence also in this \ow-px region, with an average 
relative factor of only 0.85±0.05 where factors of 0.5 - 0.6 should be expected, see also the 
discussion in Sect. [TT] below. 

Notwithstanding the apparent problems with the K"*" measurements, the interpolated data 
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at 24 GeV/c beam momentum may be integrated in order to obtain pr integrated invariant cross 
sections, mean transverse momenta and total kaon multiplicities. The resulting pt integrated 
invariant xp distributions and mean transverse momenta are presented in Fig.[68]in comparison 
to the NA49 results. 
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Figure 68: F(24 GeV/c) and {px) as a function of xp compared to the NA49 (dashed lines) and 
^/s = 3 GeV (dotted lines) results; a) and b) for K+, and c) and d) for 



For the px integrated distributions, Fig. [68k. the approach of the lower energy data to the 
NA49 results for K+, to within 10% atXF = 0.3, confirms the statements made above concerning 
s-dependence. Also the behaviour of (pr) for K+ and K~, Fig. [68b and d, raises questions, in 
particular if compared to the results at y/s = 3 GeV also shown in Fig. [68b. 

The total integrated kaon yields at y/s = 6.84 GeV, as they result from the data interpo- 
lation, are 



{n^+) = 0.107 

(24) 

(n^-) = 0.0262. 

This is 10% above and 60% below the values fitted by Rossi et al. [T| for K+ and K~, 
respectively. A more detailed discussion of total yields is given in Sects . [TTl and [T3] below. 

10.3 Data at Serpukhov energy 

In the range of 30 to 70 GeV/c beam momentum accessible at the Serpukhov accel- 
erator only a single double differential measurement of kaons is available at 70 GeV/c [fTTIl . 
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This measurement has been performed at a constant lab angle of 160 mrad which corresponds 
approximately to Xi? = and in a transverse momentum range of pt > 0.46 GeV/c. In conse- 
quence there is no possibility to establish a reliable pt integrated yield at xp = 0, not to speak 
of the total production cross section. A comparison with the NA49 data has been performed 
taking account of the dependence of Xi? on pt shown in Table |7] together with the invariant 
cross sections 
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0.343 



Table 7: Relation between px and xp for [[TTII 



The data of NA49 have been interpolated to the respective pxlxp combinations, see 
Tablelll in order to obtain the ratios of invariant cross sections presented in Tableland Fig.[69l 
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Figure 69: Cross section ratio R of flHI to NA49 for K+ and K~ 

If one interpolates the cross section ratios as shown in Fig. [69] and if one takes the 
courage to extrapolate these curves down to = GeV/c as also shown in this Figure one may 
obtain the invariant cross sections at i/i = 1 1.5 GeV, xp = 0, from the ones of the interpolated 
NA49 data and integrate over pt- This yields the invariant cross sections (Eq.[T5l) 
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F{K+,xf = 0) = 0.549 
F{K~,xf = 0) = 0.322. 



(25) 



These values are plotted in Fig. |70] together with the cross sections determined in 
Sects. 110.11 and 110.21 above, with the NA49 data and the lower range of ISR energies (see 
Sect. I10.4l below) in order to get a first view on s-dependence. 




^ [GeV] 

Figure 70: F(xf = 0) for K+ and as a function of ^/s for 6 different energies 

It is interesting to compare this dependence to the pr integrated cross sections F{xf) 
and da/dy{y) available from Zabrodin et al. [66] for only, at the lower Serpukhov energy 
of 32 GeV/c beam momentum or y/s = 7.85 GeV. They give F{xf = 0) ~ 0.6 and dn/dy{y = 
0) = 0.066. These values are about 30% higher than the ones obtained by NA49 where on the 
other hand a decrease by 2.8 would be expected from the s-dependence, Fig.|70l It is therefore 
concluded that (contrary to the pion and baryon cross section provided by ['66'|) their extracted 
yields are flawed, especially as the integration over their rapidity distribution gives a total 
yield of 0.21, about 60% higher than the value from NA49 at = 17.2 GeV/c. 

10.4 Data at ISR energy 

The ISR data on kaon production may be separated into three regions of xf- A first 
region at = is covered by [|2Tl[22l, the region from xf = 0.08 to 0.49 by EOl . and finally 
the data of lT5]-[T9ll reach from Xi? ~ 0.2 to 0.7. For the purpose of the present work these data 
are exploited in a phase space region (see Fig. [T]) in pr up to 1.9 GeV/c and in up to 0.6 in 
order to allow for a comparison to the NA49 data with reasonably small extrapolations. This 
makes available a substantial set of 383 points for K+ and 335 points for K~. 

10.4.1 The central region, l{I7\\22\I 

The data of Alper et al. [21] and Guettler et al. [l22ll follow each other with several 
years difference. The later data Il22l are extending (and superseding) the earlier work ll2Tll at 
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low pt, in a range from 0.123 to 0.280 GeV/c. They feature statistical errors in the 5 to 10% 
range, exceptionally small for ISR standards, and are probably the best controlled data as far as 
normalization and internal consistency are concerned. The combined data sets are presented in 
Fig. Ell at the five standard ISR energies from y/s = 23 to y/s = 63 GeV. 
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Figure 7 1 : Invariant cross sections at x^? = from 021112211 as functions of pr at five ISR energies, 
a) to e) for K+, f) to j) for K~. The data interpolations are superimposed on the data points 



In order to eliminate some of the larger fluctuations in the Alper et al. data [|2ni . a 
multistep eyeball interpolation imposing smoothness both in px and in the s-dependence has 
been performed, again (see Sect. 16.21) avoiding any kind of arithmetic fitting. The resulting 
Pt dependences are superimposed on the data in Fig. |7T1 The distributions of the differences 
between data points and interpolation, normalized to the statistical errors, are shown in Fig. |72] 
separately for llITIl and [|22]|. 




Figure 72: Normalized differences between data and interpolation for a) [22J b) ETTl 
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The Gaussian fit to tlie differences shows an upwards shift of 0.18 or about 1% and 
an rms of 0.6 for the data of Guettler et al. [22], indicating a certain overestimation of their 
statistical uncertainties. For Alper et al. [21] the upwards shift is equivalent to 3-5%, with an 
rms compatible with unity. The accumulation of entries at A/a in the region +0.5 and +1.5 
corresponds to the data points in the low pr region of the Alper et al. data [l2T1l visible in 
Fig. r/Tb). e), g) and j). These points are in clear disagreement with the later precision data. 
Other points of ETl . deviating far below the interpolation (region of A/a < -2) and partially 
even falling below the NA49 data, are visible notably in Fig.rTTb). c), and f). This demonstrates 
again a certain instability in the absolute normalization of the earlier ISR data also visible in 
Sects. [T0^.4.2l and [l0.4.3l and discussed for protons in [2]. 

In dividing the interpolation of lEHISll by the one for NA49 (Sect. 16.21) one may define 
the ratios 

RUx. = 0, PT, V^) = f.M9^,^ = o,pr,n.2GoY) ^'^^ 

shown in Fig.|73]as a function of pr for the five ISR energies. For comparison, also the corre- 
sponding ratios for the inter/extrapolation of the Serpukhov ifTTI and PS [SJ data are included. 




Figure 73: Rmtipr, y/s) as a function of px including data from lISlfTTI 



A remarkable picture emerges. Evidently there is a very strong change of the s- 
dependence with pr, with three clearly distinguished regions of pr- A first, low pt region 
extends up to pr ~ 0.6 GeV/c. The strong increase with s already stressed in lf22l as "rising 
plateau" is completely concentrated in this limited area. A second region at 0.6 < pr < 1 GeV/c 
shows in contrast a rather small s-dependence, limited here to a relative increase of only 10% 
(20%) for K+ and K^, respectively, over the complete range from y/s = 17 to 63 GeV. A third 
region at > 1.2 GeV/c shows again a strong s-dependence with increasing px up to factors 
of 2 (3.6) for K+ and K", respectively, over the before-mentioned range of ^/s. These fea- 
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tures reflect in the inverse slope parameters of the niT distributions (see Sect. 16.41) presented in 
Fig. 171 
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Figure 74: Inverse slopes of the interpolated rriT distributions as function of — itik for the 
5 ISR energies and for ( panels a) to e) ) and for ( panels f) to j) ) 



In plotting the extracted inverse slopes at fixed pr as a function of y/s. Fig. |75l and 
extending the s-range to Serpukhov and PS energies, the strong evolution of this "hadronic 
temperature" both with px and with and thereby the sense (or, rather, non-sense) of thinking 
in terms of a fixed "temperature" in soft hadronic production, becomes evident, see also Sect. 
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Figure 75: Inverse slopes at fixed px as a function of ^/s for a) K+ and b) K 
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below. Only (by accident) in the region of Serpukhov energies there is a concentration of inverse 
slopes in a small interval around 180 MeV for K+ and 160 MeV for K~. 



10.4.2 The intermediate x p region, / [20l/ 



In the following Sects. 110.4.21 and ll0.4.3l the NA49 data are compared to ISR results at 
Xi? 7^ 0. In this comparison, in addition to the ratio 



R{xf,Pt, Vs) 



f^^^'^{xF,PT, 17.2 GeV) 



(27) 



the ratio Rmtixp = 0,pr, a/s) (Eq. |26l) which describes the s-dependence at = is used 
in order to make a prediction at all ^ 0. As shown in these two sections, the ISR data in 
forward direction are well described by the NA49 results multiplied by -Rint(a;F = 0, pr, ^/s) at 
all Xi?. This non-trivial result shows that the s-dependence has no major change with x^?. 

The data of Capiluppi et al. [|20ll cover the ranges of 0.08 < xp < 0.5 and 0.2 < 
Pt < 1.5 GeV/c, both as a function of pr for fixed xp and as a function of xf for fixed pp. 
An overview over the pr dependence is presented in Fig. |76]for and in Fig. |77]for K^. In 
both Figures the invariant cross sections for the three xp values 0.08, 0.16 and 0.32 are plot- 
ted separately for the three \/s values of 31, 45 and 63 GeV. Also shown are the NA49 cross 
sections at these xp values and, in addition, the evolution of the cross sections at Xi? = with 
respect to NA49, see Fig. |73l for these y/s values, as a function of pr- 
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Figure 76: Invariant K+ cross sections [|20l as a function of pp at fixed xp for a) y/s = 31 GeV, 
b) ^/s = 45 GeV and c) y/s = 53 GeV in comparison to the NA49 data (dashed lines) and to the 
evolution withp^ and y/s as measured at xp = (solid lines). The results at xj? = 0.16 and 0.32 
are multiplied by 0.1 and 0.01, respectively, for better separation 



Evidently the s-dependence at x^? = is also describing the evolution in the x p region up 
to 0.3 within the sizeable statistical errors of typically 15-30%, with some exceptions notably 
for K"*" at Xi? = 0.32. The additional systematic uncertainties which can reach the same size as 
the statistical fluctuations have to be taken into account here. 
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Figure 77: Invariant cross sections f20\ as a function of pr at fixed xp for a) y/s = 31 GeV, 
b) y/s = 45 GeV and c) y/s = 53 GeV in comparison to the NA49 data (dashed lines) and to the 
evolution withpx and y/s as measured at = (solid lines). The results at = 0.16 and 0.32 
are multiplied by 0.1 and 0.01, respectively, for better separation 



A similar picture emerges for the data sets obtained at fixed px as a function of Xi;'. Here 
the ratio to the NA49 data, averaged over the xp ranges of [|20ll from 0.1 to 0.4, is presented in 
Figs.|78]and|79]as a function of ^/s, for the four values 0.21, 0.42, 0.82 and 1.27 GeV/c. 




Figure 78: Cross section ratio R with respect to the NA49 data at a) px = 0.21 GeV/c, b) 
Pt = 0.42 GeV/c, c) pr = 0.84 GeV/c, d) px = 1.27 GeV/c, averaged over the xp ranges of ll20ll . 
as a function of y/s, for K+. Superimposed is the s-dependence measured atxp = 



In both Figures the s-dependence extracted at these px values for x^ = 0, Fig. |73l is 
shown as the full line. Again the data follow this s-dependence within their statistical uncer- 
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Figure 79: Cross section ratio R with respect to the NA49 data at a) pt = 0.21 GeV/c, b) 
Pt = 0.44 GeV/c, c) pt = 0.83 GeV/c, d) px = 1.25 GeV/c, averaged over the xp ranges of [|20ll . 
as a function of a/s, for K^. Superimposed is the s-dependence measured at xp = 



tainty. 

In conclusion it may be stated that the data of EOl in the intermediate xp range from 
0.08 to about 0.4 are reasonably well described by the NA49 data supplemented with the s- 
dependence extracted at xp = 0. 



10.4.3 The forward data ofAlbrow et al. / [TTHJPl/ 



The CHLM collaboration has produced rich data sets for pions [T| and protons f2\, in the 
latter case with far more than thousand cross section values. For kaons, however, the situation 
is less favourable. In fact only less than 100 data points for each charge fall into the xp region 
below 0.5 usable for comparison purposes. On the other hand there is good overlap with the 
data [20] thus allowing for meaningful cross checks although the distributions of the statistical 
errors, Fig.[80l show wide spreads around mean values of about 15%. 
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Figure 80: Distributions of the statistical errors for the data of [fT5] - IT9l for a) K+ and b) K 
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For K~ three data sets are available. A first set [15] covers, at fixed cms angle, apr range 
from 0.16 to 0.7 GeV/c in an Xi;' window from 0.12 to 0.5, the upper cut-off being imposed here 
by the range of the NA49 data. The relation between xp and px is given by p-r = 1.33 xp. The 
ratio to the NA49 data is shown in Fig. [8T]as a function of pp separately for the three available 
cms energies averaged over x^. 
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Figure 8 1 : Ratio R between data from [fT5l and NA49 as a function of px for a) ^/s = 31 GeV, 
b) ^/s = 45 GeV and c) ^/s = 53 GeV. The pp dependence atxp = 0. Fig. |73l are included as full 
lines. The broken lines in panels b) and c) indicate the result of a 15% downwards normalization 
error 

The corresponding pp dependences at x^? = are given as full lines in Fig. [8TJ As in 
the case of the data from GOll . see Sect. 110.4.21 the general downward trend of the ratio with 
increasing pp is well described by the ratios at xp = 0, although for y/s = 45 and 53 GeV the 
data fall below (in contrast to [20]) by about 15%. This order of magnitude is definitely within 
the normalization errors typical of ISR data, as discussed in some detail in [2]. The averaged 
ratios over 0.2 < pp < 0.7 GeV/c and 0.15 < x^? < 0.5 are given in Fig. [82] as a function of 
y/s, indicating again the trend at Xi? = as the full line and the reduced ratio corresponding to 
a 15% normalization error as the broken line. 
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Figure 82: Ratio (R) between |[T5l and NA49 averaged over the intervals 0.2 < pp < 0.7 GeV/c 
and 0.15 < xp < 0.5 as a function of ^/s. Full line: behaviour at Xi;' = 0, broken line: 15% 
normalization error 

The second data set for [16] is obtained at the fixed xp value of 0.19 in the pp 
range from 0.14 to 0.92 GeV/c, at = 53 GeV. The ratio to the NA49 data, Fig.[83l shows a 
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structure which is very probably of systematic origin and equal for K~ and K+. At pr below 
about 0.5 GeV/c the ratios are compatible with no s-dependence from 17.2 to 53 GeV, whereas 
for pt above 0.6 GeV/c the values are compatible with the pt dependence observed dXxp = 0, 
full line in Fig. [831 
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Figure 83: Ratio R between fT6^1 and NA49 for K~ at xp = 0.19 as a function of px- The full 
line gives the behaviour at Xi? = 

The third and last data set available for [ 17] at ^/s = 45 GeV covers the high x p range 
above 0.5 for the three pr values 0.4, 0.55 and 0.75 GeV/c. Here the comparison is extended 
up to Xi? = 0.59 using a slight extrapolation of the NA49 K data. The ratio between [fTTl and 
NA49 has been averaged over this xp window and is shown in Fig. [84] as a function of pt- 
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Figure 84: Ratio (R) between ffTTl and NA49 as a function of px averaged over the xp 
region from 0.5 to 0.59. The full line gives the behaviour at xp = 0, the broken line shows the 
consequence of a 20% normalization error 

Remarkably, the ratio turns out to be at or slightly below unity. This would again indicate 
no s-dependence between y/s = 17.2 and 45 GeV, as compared to the behaviour at xp = 
characterized by the full line in Fig. [84] This highly improbable case could be explained by 
a 20% normalization uncertainty, see the broken line in Fig. [841 After all it should be kept in 
mind that the s-dependence in the intermediate px range where most of the comparison data are 
found is rather small and 10% effects may make all the difference in interpretation. 

For K+ four sets of data [[T6l - fT9l . with only partial overlap with the K~ results, are avail- 
able. A first set [[HI is obtained at fixed xf = 0.19 in thepp interval 0.14 < pr < 0.92 GeV/c at 
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^/s = 53 GeV. As for K^, the ratio between lfT6l and NA49, Fig. [85l shows a structure indicat- 
ing systematic problems, with no s-dependence below pr = 0.6 GeV/c followed by an increase 
of about 50% at pt ~ 0.75 GeV/c bracketing the behaviour at xp = shown as a full line in 
Fig. [851 
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Figure 85: Ratio R between (W\ and NA49 for K+ at Xi^^ = 0.19 as a function of px- The full 
line gives the behaviour at Xi;' = 



A second data set (TW\ has been obtained at constant pt = O.S GeV/c in a range of xp 
from 0.23 to 0.8 at ^/s = 45 GeV. The xp range has been cut at the upper limit of 0.6 for 
comparison purposes with the (partially extrapolated) NA49 data. The ratio between [18] and 
NA49 is shown in Fig. [86] as a function of xp. 
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Figure 86: Ratio R between [18] and NA49 for K+ as a function of xp at pp = O.S GeV/c and 
y/s = 45 GeV. The dashed line shows the mean ratio of the measured points and the full line 
shows the ratio atxp = 



The ratio averaged over the seven xp values shown is 1.43 as compared to 1.09 as mea- 
sured atxp = 0. This large value is in internal disagreement with the other CHLIVI measurements 
discussed here. 

The third data set |[T9l contains data at fixed xpinpp ranges from 0.4 to about 1 .7 GeV/c, 
at = 31, 45 and 53 GeV At = 53 GeV the four xp values of 0.3, 0.4, 0.5 and 0.6 are 
available, whereas ^/s = 31 and 45 GeV are limited to Xi;' = 0.6 only, a bit uncomfortable with 
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respect to the range of the NA49 data. Therefore only the data at ^/s = 53 GeV are compared 
here, as shown in Fig.[87l 
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Figure 87: Ratio R between [fT9l and NA49 as a function of px at a.) xp = 0.3, h) xp = 0.4, c) 
xp = 0.5, d) xp = 0.6. The pp dependence at = is shown as the full line in each panel 



Within the large error margins of (19] the ratios are compatible with flat pp dependences 
and mean values of 1.52, 1.25, 1.16 and 1.05 GeV/c for xp = 0.3, 0.4, 0.5 and 0.6, respectively. 
Except for xp = 0.3, they are also compatible with the dependence at xp = 0, shown as solid 
lines in Fig.[87J The large mean ratio at Xi;' = 0.3 is in contradiction with [|20] - l22l . 

Finally, the K+ data of [flTl at ^/s = 45 GeV cover the range inxp from 0.5 to 0.61 1 for 
Pt from 0.35 to 0.93 GeV/c. The ratio between (Tl] and NA49 is shown in Fig. [881 averaged 
over the relatively small xp window, as a function of pp. 

^ 1.8 
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Figure 88: Ratio {R) between [fTTl and NA49, averaged over xp, as a function of pp- The pp 
dependence at = is given as the full line 
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Again the dependence is compatible with a constant ratio at 1.35, but incompatible for 
Pt > 0.5 GeV/c with the behaviour at = IfTSUTOl . 

In conclusion of this sub-chapter concerning the CHLM data, a certain frustration over 
the apparent sizeable systematic effects contained in these data should be admitted. This pre- 
cludes a definite statement about the xp and pr dependences in the medium to forward region 
of longitudinal momentum. This is the more regrettable as no experiments are in view to pro- 
duce new, more precise data, especially not at the high energy colliders including of course the 
LHC. Nevertheless it may be stated that the observed patterns are compatible, taking all data 
of Sects. 110.4.21 and 110.4.3] together, with the behaviour observed at xp = within the given 
statistical errors, allowing also for the known systematic uncertainties (see also |l2l). 



10.4.4 Extrapolation ofSPS andlSR data to ^/s = 200 GeV 

In view of the scrutiny of kaon production at cms energies above the ISR at RHIC and 
the p-i-p colliders in the following sub- sections, and in view of the evident problems encountered 
with these higher energy data, it seems indicated to perform an extrapolation of the combined 
SPS and ISR data at least to RHIC energy, that is, to y/s = 200 GeV. This attempt looks feasible 
given the dense coverage of the ^/s scale between 17 and 63 GeV and the smooth behaviour of 
the cross sections as a function of y/s. This is evident from Figs. [89l and l90l where the ratios of 
kaon densities per inelastic event at = 



R' 



if {Xf,Pt)/ (^inel) 
ifixF,PT)/crmel 



ISR 



NA49 



R- 



NA49 
inel 



a. 



ISR 

inel 



(28) 



are shown as a function of y/s for fixed values of px including an extrapolation to ^/s = 200 GeV 
for K+ and K^, respectively. This extrapolation is extending the eyeball fits to the lower energy 
data presented in Sect. 110.4. ll without using arithmetic formulations. 
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Figure 89: R' as a function of for K+ at a) = 0. 1 ^ 0.8 GeV/c and b) pr = 0.9 ^ 2.0 GeV/c. 
The values of y/s are indicated with dotted lines. The NA49 point is marked with circle 



The extrapolation is facilitated by the fact that over most of the pt range covered the 
dependence on y/s is approximately linear in the double-logarithmic plots of Figs. [891 and [901 
which means a power-law behaviour of R' as a function of y/s. Noteable exceptions from this 
simple behaviour are visible both at low pt < 0.5 GeV/c and at high p^- > 1.5 GeV/c. In both 
regions the energy dependence flattens out with increasing y/s. 
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In this context it is also interesting to look at the ^/s dependence of R' towards lower 
energies which is shown in Fig. |9T]for a few px values down to = 3 GeV for K+ (Sect. 110.11) 
and 6.8 GeV for (Sect.flOj). 
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Figure 9 1 : i?' as a function of ^/s for a) K+ and b) K at pt = 0, 0.5, 1 , 1 .5, 2 GeV/c. The values 
of y/s are indicated with dotted lines 



As R' must vanish at production threshold, the wide spread of the dependences for 
constant px below SPS energies indicates a corresponding and very characteristic spread of 
kaon production thresholds with transverse momentum. This spread is charge dependent and 
reaches from ^/s about 2.5 to 10 GeV for K+ and from about 5 to 10 GeV for K^, for px from 
to 2 GeV/c. The lower effective threshold for K+ is following from the prevailing associate 
kaon-hyperon decays of non- strange baryonic resonances at low ^/s whereas K~ can only stem 
from heavy strange hyperons or heavy meson decay corresponding to a higher overall mass 
scale of the resonances involved. 

The distributions of the invariant cross sections at Xi;' = as a function of pr at 
^/s = 200 GeV, as they are resulting from the extrapolation shown in Figs. [89] and |90l are 
presented in Fig. [92] for K+ and K^, respectively. 
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Figure 92: /{xp = 0, pt) as a function of px extrapolated to y/s = 200 GeV for a) K+, b) K 

10.4.5 Pt integrated kaon yields in the ISR and RHIC energy range 

In view of the statements concerning the ISR data made above, it might seem rather 
daring to attempt the integration of the available double-differential cross sections into px in- 
tegrated and even total kaon yields. Several facts encourage, nevertheless, a new attempt based 
on purely experimental considerations: 

- The NA49 data offer a relatively precise starting point in the neighbourhood of the low- 
est ISR energy. 

- Following the discussion in the preceding sections, the evolution of the double differ- 
ential cross sections from SPS to ISR energies may be considered as experimentally 
established within error limits of about 10-30%, depending on the xp range under study 

- It is interesting to compare the integrated yields of charged kaons to the ones of K^, the 
latter ones being rather precisely determined well into the ISR energy range by bubble 
chamber experiments, see Sect. [TT] below. 

- The extrapolation to RHIC energy described above will permit a comparison of the pt 
integrated results and an estimation of the total kaon yields at ^/s = 200 GeV. 

The following approach has been followed. The detailed dependence of the invariant 
cross sections on pr and y/s established at xp = 0, Sect. 110.4.11 as characterised by the two- 
dimensional set of factors Ri^i relative to the NA49 data. Fig. |73l has been extended to the 
full range of xp. This is motivated by the comparisons with all available data discussed in 
Sects. [TO.4.21 and riO.4.31 above. This means that the invariant cross sections at each energy are 
obtained from the double differential NA49 data as follows: 



A set of cross sections at each ISR energy covering the major part of the available phase 
space is thus established allowing the extraction of pr integrated yields and total kaon multiplic- 
ities. It should be stressed here that this approach avoids the use of arithmetic formulations [jll 
which would introduce systematic uncertainties beyond the statistical and systematic fluctua- 
tions of the data. 

The resulting pp integrated quantities F, dn/dxp and {pr), see Sect. 19. II for the defi- 
nitions and the results from NA49, are presented in Table [8] for K+ and in Table |9] for K , as 
functions oixp for the five ISR energies and the extrapolation to = 200 GeV. 



f{xp,pT, y/s) = f{xp,pT, y/s = 17.2 GeV) R{xp = 0,Pt, y/s) 



(29) 
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Table 8: pr integrated quantities F, dn/dxp and (pt) as functions of xp for K+ at the five ISR 
energies and the extrapolation to = 200 GeV 
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Table 9: pr integrated quantities F, dn/dxp and {pt) as functions ofxp for K at the five ISR 
energies and the extrapolation to y/s = 200 GeV 
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Figure 93: F, dn/dxp and {pt) as functions of xp for K+ for y/s = 23, 63 and 200 GeV 



Corresponding plots of these quantities are shown in Figs. [93] and |94] for the three en- 
ergies 23, 63 and 200 GeV and for K+ and K^, respectively. Salient features of these results 
are the relatively slow and smooth increase of F with energy in comparison to the fast increase 
of dn/dxp at low xp which is practically proportional to ^i, see Eq. [15] above, and the quasi- 
invariance of mean pt with energy. The latter feature is explained by the increase of the invariant 
cross sections both at low px and at high px, compensating each other for the mean value, and 
the relatively small s-dependence in the intermediate px region. 




Figure 94: F, dn/dxp and {px) as functions of xp for K for = 23, 63 and 200 GeV 
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10.4.6 Total kaon yields 

Integration over a;^ of the dn/dxp distributions results in the total multiplicities given 
in Table [lO] together with the mean kaon yields and the total K+/K ratios. 
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Table 10: Total kaon multiplicities, mean charged yields and K+/K ratio at ISR energies and 
extrapolation to ^/s = 200 GeV 

These multiplicities are plotted as functions of in Fig.|95k. They are supplemented in 
Figs.l95b and [95b by the quantities Fixp = y/s) dud dn/ dxpixp = 0, respectively. 
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Figure 95: a) {uk), b) F{xf = 0) and c) dn/dxpixp = 0) as functions of y/s for and K 

It is difficult to define an error estimation for these quantities. In fact all values with the 
exception of the bubble chamber experiment [63] and the NA49 data for which the systematic 
errors are within a bracket of 2-12%, see Table [T] above, have been obtained using rather impor- 
tant inter- and extrapolations. It would therefore be advisable when performing comparisons or 
predictions, in particular in connection with heavy ion interactions, to allow for error margins 
of at least 20% both at energies below and above the SPS range. 

10.5 Data at RHIC 

Only rather limited experimental information is available to date from RHIC as far as 
double differential inclusive cross sections for identified kaons in p+p interactions are con- 
cerned. Data on central production come from STAR Ii23] - [25l using different identification 
methods and from PHENIX [26], both at y/s = 200 GeV, as well as preliminary data from 
PHENIX [27] at = 62.4 GeV The BRAHMS experiment has shown data a.t ^/s = 200 GeV 
with rapidities ranging from to 3.3. In the present comparison two sets of central BRAHMS 
data [l28l|29l and the most forward data at rapidity 2.95 and 3.3 fl30l will be addressed. 
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The invariant cross sections at y = and y/s = 200 GeV, Fig. [96l form a wide band 
within a margin of about a factor of 1.5-2 in the pr range from a lower limit at 0.25 GeV/c for 
STAR and about 0.4 GeV/c for PHENIX and BRAHMS up to the upper hmit at about 2 GeV/c 
usable for the direct confrontation with the lower energy data in this paper. 




Figure 96: Invariant kaon cross sections from RHIC at central rapidity and y/s = 200 GeV [|23] - 
|26l|28l|29l as a function of pr for a) K+ and b) K~. The corresponding data from NA49 (dotted 
line), the interpolated data from the ISR at y/s = 63 GeV (dashed line) and the extrapolated ISR 
data at y/s = 200 GeV (full line) are also shown 

The reason for this large internal variation of results which goes beyond any other data 
sets discussed in the preceding sections, has to remain open for the time being. It is however 
clear that, if compared to the NA49 data, to the ISR data at y/s = 63 GeV and to the extrapolated 
ISR data. Sect. 110.4.41 also shown in Fig. [961 there is an evident change in shape of the p^ 
distributions in particular compared to the 200 GeV extrapolation. In the lower region, below 
about 1 GeV/c, all RHIC data approach or cut below the ISR extrapolation, whereas towards 
high pt a rather constant, large offset of factors 2 to 3 is visible. 

This is quantified in the ratio plots shown in Fig. |97] Here the ratios R' of particle 
densities per inelastic event. 



r/ \ I v/i=200GeV 

f{ XF,PT)/K,l 

f{xF,PT)/crl^i 



^ - 77 N / NA49 



are presented in order to take out the increase of the inelastic cross sections with energy for the 
PHENIX data [26] and the extrapolated ISR data (Sect. 110.4.4]) 

In Fig. [971 two basic features of this data comparison are clearly visible for the PHENIX 
data: both for K+ and for there is at > 1 GeV/c a nearly constant factor of 1.5-1.8 with 
respect to the extrapolated ISR data, whereas for pr < I GeV/c the data sets approach each 
other rapidly to become equal at the lower pt cut-off of the RHIC data. 
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Figure 97: Ratios R', Eq.[30l as functions ofpr, for the data from PHENIX ^ (full line) and 
the extrapolated ISR data at = 200 GeV (dashed line). Panel a) for K+, b) for 

At this point it might be useful to look at the central kaon data from PHENIX at 
^/s = 62 A GeV |l27l that is, in the immediate neighbourhood of the ISR data lIHllIIll at 
y/s = 63 GeV. Here, the PHENIX experiment gives a px distribution with the same lower cut-off 
as at 200 GeV, at pr = 0.45 GeV/c as shown in Fig. 1981 
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Figure 98: Invariant kaon cross sections at ^/s = 62.4 GeV from PHENIX [|27ll . for a) K+ and 
b) K^, as a function of pt- The original data points (full points) and the data divided by vr (open 
points) are indicated. The data from the ISR [I2T1I22]| at ^/s = 63 GeV are also shown, together 
with the NA49 data, as the full and dashed lines, respectively 

A glance at the RHIC data points in Fig. |98] shows that they are far above the ISR data by 
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a factor of 3-5 varying with px, and even well above the PHENIX data at ^/s = 200 GeV [|26l . 
As the same large factors apply for pions and baryons, it has been concluded that a factor of 
I/tt has probably been dropped in the cross sections given in [|271l . a fact that is not uncommon 
in the definition of rapidity densities. Tentatively applying this factor, the cross sections move 
down to the lower data points shown in Fig. |98l In direct comparison to the ISR data one may 
define the ratio 



e(pr, xf = 0) 



shown in Fig. 
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Figure 99: Cross section ratio ei^px) as a function of px for a) K+ and b) K^. The full lines 
give a representation of e as a constant offset at pt > 1 GeV/c combined with an efficiency loss 
towards lower px 

Both for and for K+ two features are emerging from this plot. Below px ^ I GeV/c 
there is a sharp drop of e reaching values below 1, corresponding to cross sections below the 
ISR data at the lower px cut-off at 0.45 GeV/c. This looks like an apparative loss of efficiency 
for kaon detection towards low px- At px > ^ GeV/c on the other hand there is an offset which 
is approximately px independent at a value of about 1.3 for and 1.45 for K+. Tentatively 
regarding the ISR data as a reference this may be translated into a correction factor to be applied 
to the PHENIX data as a function of px indicated by the full lines in Fig. |99] which combine 
a constant offset determined at px > I GeV/c, with an efficiency drop towards lower px- The 
latter effect, if of apparative origin, might be expected to hold for all reactions and all interaction 
energies studied by this experiment, in particular also for the measurements at y/s = 200 GeV 
both for elementary and nuclear collisions. The overall offset, on the other hand, could well 
depend on different experimental constraints as for example vertex distributions and/or trigger 
efficiency, and thereby be s and reaction dependent. In particular the trigger conditions are 
largely different for elementary and nuclear reactions. One critical factor in comparing p+p 
interactions between ISR and RHIC experiments is given by the fraction of inelastic events 
picked up by the trigger arrangements, with trigger efficiencies approaching 100% at the ISR 
as compared to typically 60-70% at RHIC which favours small impact parameters and thereby 
will tend to enhance strangeness yields. 

Coming back now to the situation at y/s = 200 GeV, Fig. |97l one may try to apply the 
correction factor e{px, xp = 0), Eq.[3Tl as determined from the PHENIX data at y/s = 62.4 GeV, 
to the higher energy data, allowing only for an additional constant overall factor corresponding 
to a variation of the offset term. As shown in Fig. 1 1001 an additional factor of 1.3 applied to 
e{px,XF = 0) both for and brings the RHIC data sets into close agreement, within a 
10% margin, with the extrapolated ISR data. 
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Figure 100: Ratios R' as functions of px, for the data from PHENIX [26] divided by 1.3-e (full 
line) and the extrapolated ISR data at ^/s = 200 GeV (dashed line). Panel a) for K+, b) for K~ 



This admittedly rather daring procedure might nevertheless bring some consistency into 
an experimental situation which otherwise would appear distressingly incoherent within large 
factors. 

As far as the STAR results are concerned, they seem to indicate a similar combination 
of droop at low px and a constant, very large overall offset. In view of the sizeable internal 
inconsistencies between the different publications from this experiment, a comparable study 
has however not been tried here. 

The central BRAHMS data follow the PHENIX cross sections rather closely for K+ 
down to a lower cut-off in px at 0.55 GeV/c [|29ll . Concerning the K+ and K data shown 
in ll28l there is however a rather dramatic and unphysical drop at the given lower px limit at 
0.375 GeV/c indicating an efficiency loss very similar to the one observed for PHENIX. In [|28l 
the data for both charges fall below the PHENIX values by about 20% in the overlapping px 
region. 

The final data sample addressed in this comparison concerns the forward measurements 
from BRAHMS at rapidities 2.95 and 3.3 and y/s = 200 GeV [|30l. As shown in the xp/pr 
correlation plot of Fig. IIOII these data start from a lower limit at about 0.7 GeV/c in px and 
correspond to an x^? range between about 0.1 to 0.3 for the combined two rapidity values. 




"0 0.1 0.2 0.3 0.4 

Figure 101: Correlation between xp and px for the two rapidities 2.95 and 3.3 of ll30l at 
^/s = 200 GeV/c 
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This coverage is comparable to the intermediate data of Capiluppi et al. [|20| at the ISR, 
see Sect. I10.4.2I These results have been shown, Figs. |76] to |79l to be compatible with the 
application of the s-dependence observed at a; = to the NA49 data in the corresponding xp 
and pt ranges. It is therefore interesting to confront the forward BRAHMS data both with the 
NA49 data and with the extrapolation of the ISR results to a/s = 200 GeV, Sect. I10.4.4[ as 
shown in Fig. I102[ 




Figure 102: Invariant cross sections as a function of px for rapidity 2.95 and 3.3 for K"*" (panels 
a and b) and (panels c and d), respectively. Also shown are the NA49 data (dashed lines) 
and the extrapolated ISR data (full lines) at these rapidities 



Evidently the K+ data from BRAHMS are rather close to the extrapolated ISR data for 
Pt > 1.2 GeV/c, whereas the K~ data show offsets by factors of about 0.6 aty = 2.95 and 0.5 
at y =3.3 in the same px range. Below ~ 1.2 GeV/c the BRAHMS data increase rapidly up 
to a local maximum at px ~ 0.8-0.9 GeV/c which is evidently non-physical. The sharp drop of 
the cross sections below this maximum to values even below the NA49 data indicates again the 
loss of kaon detection efficiency below ~ 1 GeV/c which seems to be common to all RHIC 
data which have been discussed in this section. These features are quantified in Fig. 11031 where 
the ratio of kaon densities per inelastic event R' (Eq. [28]) is plotted as a function of transverse 
momentum. 

In conclusion to this section it appears that the RHIC data discussed here seem to in- 
dicate not only problems with absolute normalization evident in the comparison of different 
experiments at the same energy and to extrapolations from the ISR range, but in addition a 
common drop in kaon efficiency in the approach to their lower pt cut-off, see also Sect. 110.6.11 
This cut-off is, with about 0.4 to 0.8 GeV/c, uncomfortably high with respect to an eventual 
determination of (pt) ■ The use of these data as a reference for nuclear interactions, in particular 
concerning eventual "nuclear modification" or "jet quenching" effects widely claimed by the 
RHIC community, is therefore to be seen with some concern. 
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Figure 103: Ratio R' of kaon densities per inelastic event as a function of px for K"*" (panels a 
and b) and K (panels c and d). The ratios R' corresponding to the extrapolated ISR data are 
shown as the full lines 



10.6 Data from p+p colliders 

Three experiments at the CERN p+p collider have given kaon cross sections: UA5 WJ\ - 
ETl, UA2 fTl] and UAl [J3J in the range of y/s from 200 to 900 GeV. At the Fermilab Tevatron, 
two groups, CDF [I74ll75ll and E735 [76J have produced kaon data from y/s = 300 to 1800 GeV. 
These data are generally centered at central rapidity, within a range of 1.5 to 5 units. From 
refs. [7J] and [iTGll only unnormalized yields are available. For charged kaons, the statistical 
significance is limited to a few dozen to a few hundred identified particles, whereas for a 
wide range from a few hundred up to 60k reconstructed decays is covered. Due to the isospin 
configuration of the initial state and the limited acceptance of all experiments in xp, only the 
mean charged yields, (K"*" + K~)/2 are given. As also the equality: 



K°5 = (32) 



K+ + K 

y 

is at least within the quoted errors fulfilled for this energy region, see the following Sect. [TT] for 
a more detailed argumentation, both the mean charged kaon and the data are combined in 
this section in an attempt to link the results to the lower energy regime discussed above. 

As all experiments use double-arm triggers with a limited coverage in the extreme for- 
ward direction, the trigger cross sections correspond in general not to the total inelastic cross 
section but to a fraction of the so-called "non single-diffraction" cross section. This fraction is 
quoted as 93% (E735), 95% (UA5), 96% (UAl) and 98%(UA2). Since single diffraction makes 
up about 15% of the total inelastic cross section, the experiments trigger on about 80% of ainei- If 
compared to the NA49 and ISR data including the extrapolation to 200 GeV which are obtained 
in relation to the full ciinei, a correction for the trigger losses is in principle necessary. Only the 
UA5 collaboration has estimated this correction ||69l to about -16% at y/s = 200 GeV and -12% 
at y/s = 900 GeV. In the following subsections all comparisons are carried out including the 
necessary correction to the full inelastic cross section. In addition all data are given as invariant 
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densities by dividing the invariant cross sections, if given in mb, by the inelastic cross section. 



10.6.1 Data at ^ = 200 GeV 

The UA5 experiment [711 gives cross sections for and (K"*" + K~)/2 which may be 
compared to the extrapolation from NA49 and ISR data to this energy, Sect. llO.C4l As the UA5 
data are given over their rapidity interval of ±2.5 units of rapidity in the form 1 / a^^sud^ cr / dp^ 
a transformation into invariant density 1 / {2T[pT(Tme.\)d^ a / dydpT has been performed including 
the correction for trigger losses, see above. This results in the data shown in Fig . 1 1 04b compared 
to the extrapolated NA49/ISR data. 
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Figure 104: a) Invariant kaon density from UA5 compared to the NA49/ISR data extrapolation 
(dashed line) and to the extrapolation multiplied by 1.35 (full line) b) UA5 fit to their data 
(dashed line) compared to the data extrapolation multiplied with 1.35 (full line) 

Evidently the UA5 data are on average higher than the data extrapolation (dashed line 
in Fig. ll04h ) by about 35% as shown by the full line in Fig. ll04h . This systematic difference is 
certainly compatible with the uncertainty inherent in the data extrapolation and with the ~20% 
uncertainty given for the normalization of the UA5 data. What is interesting here is that the 
shape of the extrapolated distribution after renormalization is compatible within the statistical 
errors with the UA5 data over the full range of pr from 0.07 to 2 GeV/c. 

UA5 has performed a fit to their data of the double form: 

1 da (Ae-'"^^, forpr <0.4GeV/c 

= < / \ " (33) 

a^sDdp^T U'(^) ' forpr>0.4GeV/c 

The first form at low px is necessitated by the unphysical behaviour of the second form 
through pt = 0; it is motivated by the idea of thermal behaviour at low transverse momen- 
tum. The fit parameters [|6l at = 200 GeV are A = 10.9, b = 8.2, A' = 0.60, n = 8.8 and 
Po = 1.3 GeV/c. The inverse mr slope of 0.12 GeV/c implied by the parameter b is however 
rather low and corresponds to the non-thermal behaviour observed for the ISR data, see Fig. 1741 
As stated above, the fit values have been reduced by 16% for comparison at the full inelastic 
cross section. 

The UA5 fit is compared to the renormalized data extrapolation (factor 1.35 introduced 
above) in Fig. 1104b where good agreement is visible down to pr ~ 0.5 GeV/c. The deviation 
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towards pr = GeV/c leads to a difference of about 20 MeV/c in {pt), see Sect. 110.641 below. 
As far as integrated yields are concerned, UA5 gives a rapidity density of 0.12 per inelastic 
event ||69l which corresponds to the integrated cross section F = cTmei/TT ■ dn/dy = 1.608 mb 
and compares to 1.277 mb for the NA49/ISR extrapolation. This is a 26% difference which 
agrees, taking into account the different shape of the pt distributions, with the renormalization 
shown in Fig. 1 104b . For the total yield, UA5 extrapolates to full phase space using model as- 
sumptions [l69ll . This leads to a total K^. multiplicity of 0.68 IItTII or 0.72 |[69l per inelastic event 
at y/s = 200 GeV. For the NA49/ISR data extrapolation this number is ('T'(k++k^)/2) = 0.735 per 
inelastic event. This agreement to within 5% is of course to be regarded as fortuitous in view of 
the large uncertainties involved in both attempts to estimate total yields. 

A further interesting comparison is offered by the data from STAR [l24l concerning 
invariant densities per event 1 / {27rpT) ■ d^N / dydpT- These data are shown in comparison to the 
NA49/ISR extrapolation of (K+ + K-)/2 in Fig.fTUSk. 
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Figure 105: Comparison of data from STAR at RHIC with the (K+ + K")/2 extrapolation 
at y/s = 200 GeV. Panel a) invariant rapidity densities, panel b) ratio between the two results, 
panel c) ratio between the STAR results and the UA5 data fit 

At Pt > 0.8 GeV/c there is a large offset of more than a factor of 2 between the STAR 
data and the extrapolation. This offset reduces rapidly towards lower pr until the STAR yields 
fall below the extrapolation at their lowest measured pr- The ratio between STAR and extrap- 
olation is given in Fig. 1105b . The observed behaviour with an offset factor of 2.4 and a rapid 
decrease of the ratio below ~ 1 GeV/c reproduces the features seen for charged kaons, see 
Sect. 110.51 The comparison of the STAR data with the UA5 data fit. Fig. 1105b . shows that these 
data are also in disagreement with the UA5 results obtained at the same energy. 

10.6.2 The ^ region of 540-630 GeV 

Five measurements are available in this region at i/i of 540 and 630 GeV: (K+ -i- K^)/2 
and from UA5, (K+ -i- K-)/2 from UA2 and (K+ -i- K-)/2 from E735 at = 540; (K+ -i- 
K")/2 and from UAl and K% from CDF in two different data sets at = 630 GeV. As the 
central rapidity density dN/dy changes only by 4.8% and the inelastic cross section by 1.6% 
between these two energies, the results may be compared by introducing the resulting small 
correction. In the following all results will be referred to ^/s = 540 GeV. As in Sect. 110.6.11 
invariant densities will be obtained from mb cross sections, whenever given, dividing by the 
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inelastic cross section. In addition a reduction of 14% is introduced to refer the data to the full 
inelastic cross section. 

The fit (l33l) to the UA5 data ||69ll has been chosen as an absolute reference for the sub- 
sequent data comparison. At y/s = 540 GeV the parameters are A = 7.09, b = 7.5, A' = 0.508, 
n = 7.97 and po = 1.3 GeV/c. The comparison of the UA5 data [68J with this fit is shown in 




Figure 106: UA5 data at y/s = 540 GeV. The full line corresponds to the fit described in the text 

UAl data with about 60k and 3000 charged kaons are available [|73l . These data are 
shown in comparison to the UA5 fit in Fig. 1 1071 up to pt = 2.08 GeV/c. 
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Figure 107: a) and (K+ + K")/2 data from UAl (H in comparison to the fit of the UA5 
data (full line) b) ratio /?uai = ( —) I f — ^ • The mean offset factor of 0.634 for 

V^^'/UAl/ V / UA5 fit 

is indicated in panel b) with dashed line 



For the Y^g there is a mean offset by a factor of 0.634, with an excellent reproduction of 
the shape of the UA5 fit as a function of p^. This is quantified in panel b) of Fig. [107] where the 
ratio between the UAl data and the fit is presented as a function of p^- The fluctuation around 
the mean ratio, with a standard deviation of about 7%, is compatible with the errors quoted 
by UAl. The four given data points for charged kaons, with substantially larger errors, are on 
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average higher than the K° data and fluctuate to within one standard deviation around the UA5 
fit. 

UA2 f72] has published seven data points on (K+ + K~)/2 yields which are compared 
totheUAS fit in Fig.fTOSl 
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Figure 108: a) (K+ + K-)/2 data from UA2 BH at = 540 GeV compared to the UA5 fit, b) 
ratio R\jA2 = ( — ) / ( — ) • The mean offset factor of 1.13 is indicated in panel b) 

V o-,„ei J / V / UA5 fit 

with dashed line 

As is visible from the ratio as a function of px in Fig. 1108b there is good agreement 
between the two data sets, with a mean offset of only +13% of the UA2 data with respect to the 
UA5 fit. Again the shape of the pt distribution is well reproduced. 

The CDF collaboration at the Fermilab Tevatron has published two data sets concerning 
production at y/s = 630 GeV. The first set (CDF I), with only 27 measured, yields 6 
absolutely normalized data points compared in Fig. |109| to the UA5 fit. The second set (CDF II) 
with the very large statistics of 32k is not absolutely normalized. It has been re-normalized 
to the UA5 fit at = 1 -55 GeV/c. 
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Figure 109: Comparison of K° data from CDF with the UA5 fit, a) cross sections as a function 
of Pt- CDF I data (full circles), re-normalized CDF II data (open circles); Ratios b) -Rcdfi = 
(—) I f— ) ' c) i?cDFii = I ■ The mean offset of CDF I 

V '^'"=1 / CDF I / V / UA5 fit V ^^'"=1 / CDF II / V '^'"'=' / UA5 fit 

data from UA5 fit is indicated in panel b) with dashed line 
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As quantified in the ratio plots, Fig. 1 109b and c, the absolute data CDF I fluctuate around 
the UA5 fit, with a mean offset of only a couple of percent, in the range 0.85 < pr < 2.35 GeV/c. 
The CDF II data on the other hand, which cover a larger pt range starting at 0.45 GeV/c, show 
after re-normalization to the UA5 fit large systematic deviations from the UA5 fit which increase 
sharply below ~ 1-5 GeV/c, as presented in Fig. 1109b . 

In this situation the (K+ + K^)/2 data from the E735 experiment at the Tevatron [|76l . 
although not absolutely normalized, give important information in this lower region as to the 
shape of the px distribution. Fig. II 101 
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Figure 110: Comparison of the (K+ + K )/2 data from E735 with the UA5 fit, a) data and fit as 
a function of px, b) ratio -Re735 = ( — ) / ( — ) as a function of pr- The mean offset 

V / E735 / V / UA5 fit 

factor of 0.958 is indicated in panel b) with dashed line 



The 9 data points given, after re-normalization to the UA5 fit at pt = 0.45 GeV/c, clearly 
support the shape of the UA5 fit in the region 0.2 < p^ < 1.2 GeV/c, as compared to the 
deviating CDF II data. This is quantified in Fig. II 10b with a mean deviation by a factor of 0.958 
and fluctuations which comply with the given error bars. 



10.6.3 Data at ^= 1800 GeV 

Only the CDF experiment, again with two data sets (CDF I [74] and CDF n f75\) and the 
E735 collaboration 117611 have published kaon data at the highest Tevatron energy of 1800 GeV. 
Here the fit to the CDF I data, transformed to kaon densities by dividing by the inelastic cross 
section, and corrected by -14% for the trigger losses, is used as a reference. The fit has the form 
//c^inei = C/{po+ptT with C = 5.38, n = 7.7 andpo = 1-3 GeV/c. As shown in Fig.[III]it has 
been modified at px < 0.4 GeV/c following the shape of the UA5 fit in this px region, in order 
to avoid the unphysical behaviour of this form at low px. 

The 9 data points given for the CDF I sample, corresponding to about 450 Kg, are given 
as full dots in Fig. II 1 Ih . The data points from CDF II (open circles) deviate again from the fit 
for < 1.5 GeV/c. This deviation. Fig. II 1 lb . reproduces exactly the phenomenon observed at 
y/s = 630 GeV, see Fig. 11091 :, thus indicating a systematic problem in the CDF II data analysis. 
On the other hand the re-normalized E735 data trace the CDF I fit rather well as a function of 
Px, Fig. llllb . supplementing the px scale of CDF I which is limited to px > 0.8 GeV/c, towards 
low Px. 
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Figure 111: Kaon data at ^/s = 1800 GeV; a) Full line fit to the CDF I data [741. Full circles: 
CDF I data. Open circles: CDF II data re-normalized to the CDF I fit at px = 1.55 GeV/c. 
Triangles: E735 data re-normalized to the CDF I fit. b) Ratio between the re-normalized CDF II 
data and the fit as a function of pr c) Ratio between the re-normalized E735 (K+ -i- K^)/2 data 
and the fit. The mean offset of E735 data from CDF I fit is indicated in panel c) with dashed 
line 




10.6.4 Mean transverse momenta 

Given the uncertainties and partial inconsistencies of the collider (and RHIC) data dis- 
cussed in the preceding sections, especially concerning the general lack of coverage and the 
evident systematic deviations in the low Pt region, it is not surprising to perceive large varia- 
tions in the first moments of the px distributions. Indeed, if the mean transverse momentum of 
K° or (K"*" -I- K^)/2 is plotted as a function of in the RHIC and p-i-p collider energy range, 
Fig. ll 12[ a rather disturbing overall picture emerges. 
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Figure 112: {px) in the range from 200 to 1800 GeV from different experiments 

The data which have been published in a time window from 1985 to 2008 span an 
extremely wide band of typically 0.2 GeV/c at each of the 5 available energies. For each of the 
experiments certain assumptions about the shape of the pr distributions have to be made (see 
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the preceding section for some examples) and in all cases some extrapolation either towards 
low Pt or towards high pr has to be established. 

In order to bring the evaluation of the mean transverse momentum and the data compar- 
ison on a more quantitative level, the following definitions have been used. 

The kaon density from which the mean transverse momentum is derived may be defined 
in bins of as a function ofxp and in bins of y as a function of y: 

dn , dn ,^ 

and ——. (34) 



dxpdpT dydp- 
The corresponding mean pt values are: 
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(35) 



(36) 



where E is the kaon energy and / the invariant inclusive cross section, Sect. [51 

Evidently there is a difference between the two definitions given by the energy factor in 
{Pt)xp (Eql35l). This term will enhance the contribution from low pt and reduce the contribution 
at high p-r to the mean value in Xi? as compared to the y binning. In addition at y unequal to 
the longitudinal dependence of the cross sections will couple into the mean value {pT)y as 
well as the kinematic limit in ptot which will truncate the pr distribution at small angles. The 
resulting xp and y dependences of {pt)xp and {pT)y are shown in Fig. II 13l for the case of the 
NA49 experiment at = 17.2 GeV. 




Figure 1 13: a) {pt)xp as a function of xp and b) {pT)y as a function of y at y/s = 17.2 GeV 

Clearly {pT)y is bigger than {pt)xf atxp = y = 0,at this energy by about 70 MeV/c. One 
may question the extension of (pr)y to |/ > as then a rather complex interplay of transverse and 
longitudinal dependences intervenes. The y distribution of {pT)y therefore decreases steadily 
with y whereas {pt)xp shows a characteristic increase with xp ("seagull" effect). 
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The {pt) values shown in Sect. 19. II above are defined in Feynman xp, whereas all {pr) 
values at collider energies shown in Fig. II 121 are defined in rapidity bins. In addition, the pt 
integration for the lower energy data has been established in the range < px < 2 GeV/c. 
Hence these results are not directly comparable. 

The following procedure has therefore been adopted. In the collider energy range 200- 
1800 GeV the results with doubtful cross section behaviour at low pr 112311751 are not con- 
sidered for their (pr) values. The fits to the UA5 (691 and CDF I [741 data are used at 
= 200, 540 and 1800 GeV for the determination of both {pt)xf and {pT)y For the lower 
energy data (Sect. 19.11 ) {pT)y is calculated in addition to {pt)xf for (K+ + K~)/2 including the 
data extrapolation to y/s = 200 GeV. In a first step, the integration is carried out in the range 
< Pt < 2.0 GeV/c in order to obtain for all data a comparable basis in this lower range.The 
resulting {pt)xp and (pr)?/ values are shown in Fig. II 141 as a function of y/s. 



Figure 114: {pt)xf at x^? = and {pT)y at y = as a function of -y/i. The lines are drawn to 
guide the eye 

This Figure exhibits a smooth behaviour of {pt) in the energy range from yfs =11.5 
to 63 GeV, also including the extrapolation to 200 GeV, with a variation of only 30 MeV/c for 
{Pt)xf and {pT)y There is, however, a clear offset of about 10 MeV/c for {pTjxp and 20 MeV/c 
for {pTiy between the extrapolation of the lower energy data to = 200 GeV and the trend 
of the collider data which cannot be imputed to high pt tails in this integration window. It is 
rather the different behaviour at low px, see Fig. 1104b . which can explain the difference. Given 
the general uncertainty of the collider data in the pt range below 0.5 GeV/c, the observed offset 
may still be regarded as compatible with the published errors which are on the level of 30 to 
40 MeV/c ||69ll . Another interesting feature is the rather small increase of {pt) which is only on 
the order of 50 MeV/c for {pt)xf and 70 MeV/c for {pT)y between yfs = 200 and 1800 GeV, 
always in the pt range below 2 GeV/c. 

In order to quantify the dependence of (px) on the upper integration limit in pr, this 
limit has been increased from 2 GeV/c to 6 GeV/c. For this study the published polynomial fits 
of the collider data have been used. For the lower energy data the following procedure has been 
chosen. The polynomial form 



plJP = 2 GeV/c 
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has been fitted to the high px region px > 1.5 GeV/c with po fixed at 1 .3 GeV/c. This procedure 
is possible for y/s =11.5 GeV where the data reach up to pt = 4.2 GeV/c and in the ISR energy 
range where data up to = 4 GeV/c are available. The corresponding exponents n are plotted 
in Fig. ll 15l as a function of ^/s. 
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Figure 115: Fitted exponent n as a function of ^/s. The full line is shown to guide the eye 

A consistent and smooth drop in the exponent n from about 16 at ^/s = 11.5 GeV to 8 at 
y/s = 1800 GeV is evident, describing the flattening of the pt distributions with increasing y/s. 
This allows for the interpolation ton= 14 at y/s =17.2 GeV where the NA49 data do not reach 
beyond =1.7 GeV/c. 

For the y/s range below ISR energies, the kinematic limit in px at = 2pT / yfs = 1 has 
to be taken into account. This limit influences the measured yields progressively from xt = 0.5 
upwards. This necessitates a downward correction of the polynomial fit at pt > 3 (4.5) GeV/c 
for y/s =11.5 (17.2) GeV, respectively. 

The increase of the mean px values as a function of the upper integration limit from 2 to 
6 GeV/c is given in Fig. ll 161 where the difference (pt)-{pt) up=2Gev/c is shown for y/s from 1 1.5 
to 1800 GeV both defined inxp and in y bins. 
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Figure 116: (pt)-(pt)p"p=2G'=v/c for a) {pt)xf ^^d b) {pT)y The y/s values range from 11.5 to 
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The (pt) values saturate rapidly at an upper integration limit of 2.5 to 3.5 GeV/c between 
Serpukhov and ISR energies, with a total increase of less than 3 MeV/c (10 MeV/c) for {pt)xp 
and {pT)y, respectively, in this energy range. At collider energies the saturation limit moves 
up to beyond 6 GeV/c, with very substantial increases of more than 15 MeV/c and more than 
50 MeV/c in {pt)xp and {pT)y, respectively. {pt)xp and {pT)y are shown in Fig. 11171 as a 
function of ^/s for the upper integration values from 2 to 6 GeV/c. 




The rather complex dependence of the mean pt values both on y/s and on the upper 
integration limits, in addition to the apparent systematic effects in passing from the ISR to 
collider energies, calls for some remarks: 

- A precision measurement of (pt) with an absolute error of less than 20 MeV/c in the 
region above ^/s = 100 GeV/c is still missing. This fact is mostly due to uncertainties in 
the low Pt region. 

- The sizeable difference between {pt)xf and {pT)y has to be taken into consideration 
whenever results on (px) are to be compared for different experiments and ^/s regions. 

- The large dependence of (pr) on the upper integration limit in the collider energy range, 
especially for {pT)y, is a reason for concern. It may be asked whether the definition of 
an average quantity which depends strongly on a high p-r tail more than a factor of 10 
above its value, makes any sense. 

- In fact at least part of the increase of (px) with y/s is to be imputed to the extension 
of the available transverse phase space. The kinematic limit in px is below 6 GeV/c at 
Serpukhov energy and it must be recalled that this limit is influencing the particle yields 
already at xt = 7.pTl\fs > 0.5, this means above ~ 3 GeV/c at this energy. 

- The above remarks are especially applying for the dependence of (px) on additional 
constraints, as for instance on the total hadronic multiplicity. Also in this case it might 
be advisable to separate clearly the behaviour in the lower region from the increasing 
high Pt tails. 
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11 The s-dependence of K° production and its relation to charged kaons 

A sizeable number of experiments Il39] - l44ll46]459l have addressed neutral kaon produc- 
tion from a/s = 3 GeV to y/s = 27.6 GeV. This ensures coverage from close to threshold up to 
well into the ISR energy range. Essentially all these measurements come from Bubble Cham- 
bers. This has the consequence that the total number of reconstructed is usually rather limited 
to a range between a few hundred and a few thousand. This handicap is however offset by the 
superior quality of the Bubble Chamber technique in terms of reconstruction efficiency, control 
of systematic effects and corrections, and above all a well-defined absolute normalization. It is 
in particular interesting to compare the K° to the average charged kaon yields discussed above, 
as the Eq. [32] is generally assumed to hold based on isospin symmetry [|68il although it is not 
fulfilled for instance for </> and Charm decay. 

Due to the low event statistics, double differential cross sections are not available from 
any of the experiments with the exception of [46]. The pt integrated invariant cross section F 
(see Eq. [15]) has however been given by 10 experiments between ^/s = 4.9 and ^/s = 27.6 GeV. 
These data are plotted in Fig. lllSl as a function of xp ■ 




Figure 118: pr integrated invariant cross sections F as a function of xp for various values 
of ^/s. Independent hand-interpolations at each energy are given as full lines. The dashed lines 
correspond to -F((K+ + K")/2) from NA49, Sect. [9711 Table [5] The interpolated results from 
ISR [|T5ti22]| for F((K+ + K-)/2) are also presented 
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Independent hand-interpolations at each energy have been performed in order to allow 
for the evaluation of the s-dependence at fixed values of as shown in Fig. II 191 




Ms [GeV] 

Figure 1 19: pr integrated invariant K^. cross sections F as a function of interpolated to 
fixed values of from to 0.4 using the hand-fits shown in Fig. 11181 The error bars are an 
estimation of the uncertainties of the full lines in Fig. 11181 The lines are drawn to guide the eye 

Also shown in this Figure is F((K^ + K^)/2) as derived above for the -extrapolated 
Serpukhov data at ./i = 1 1 .5 GeV, the NA49 data and the interpolated ISR data, Sect. 110.41 
Evidently there is agreement, within the experimental uncertainties, with the interpolated 
data at all xp values. This might lend some credibility to the assumptions contained in the 
evaluation of the ISR data over the full phase space in Sect. [lOl 

A similar procedure may be performed for the total integrated yields per inelastic 
event in comparison to the total mean charged kaon yields. This comparison is shown in Fig. 1 1201 




^ [GeV] \Is [GeV] 

Figure 120: Total integrated K° yields as a function of a/s Il39] - l44ll46] - [59ll (open circles). The 
total mean charged kaon yields for NA49 (full circle), ISR and extrapolation to 200 GeV (trian- 
gles) are also shown. The scale in panel b) is extended up to ^ = 200 GeV. The lines are drawn 
to guide the eye 
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using the total number of K^. per inelastic event given by the 19 experiments refs. Il39l - l44ll46] - 

m. 

Also for the total kaon yields the equality (l32l) is fulfilled within errors. Note that the 
NA49 data have a 3% error bar corresponding to the estimated systematic uncertainty, whereas 
the interpolated ISR data have been tentatively given a 10% error bar. 

12 Some remarks about contributions from resonance decay 

The evolution of the observed kaon yields with transverse momentum and interaction 
energy described in the preceding sections is characterized by rather complex patterns which 
are not easily describable by straight-forward parametrizations as they might follow from parton 
dynamics or thermal models. It seems therefore reasonable to evoke for illustration the contri- 
bution from the decay of some known resonances to the inclusive kaon cross sections. Three 
resonances, the 0(1020), the A(1520) and the charmed mesons D(1865) have been selected here 
as they give an idea about the build-up of kaon yields at low px for the two former cases, and to- 
wards high pt for the latter one. In this context it should be recalled here that most if not all final 
state hadrons are known to be created by the decay of mesonic and baryonic resonances [|35[|77ll . 
Indeed, the estimations quoted in B5[|77ll . using only a limited set of mesonic and baryonic res- 
onances, arrive at fractions of 60-80% from resonance decay for all studied final state hadrons. 
See also [34 ] for a more recent study based on two-body decays of 13 known resonances. 

12.1 0(1020) and A(1520) production and decay 

The (f) production has been measured by a number of experiments for p-i-p interactions 
in the SPS energy range 1135117811 . Results from the NA49 experiment 11791 are being used here 
to obtain the inclusive dn/dxp and cPa/dp^. distributions shown in Fig. 11211 
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Figure 121: a) a; and b) p% distributions of 

These distributions integrate to {n^) = 0.0143 per inelastic event or an inclusive cross 
section of 0.453 mb in good agreement with other measurements. 

Due to the very low Q value (32 MeV) of the (p decay into two kaons, the resulting pt 
and xp distributions are narrow compared to the inclusive kaon cross sections. This is reflected 
in the ratios Rf^^ of K~ mesons from (p decays to inclusive shown in Fig. 11221 

Evidently this contribution is very sharply peaked at small px and vanishes at xp > 0.3. 
The given percentages have to be regarded as lower limits, as production is known to be 
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accompanied by additional kaons in most if not all cases iTTSll . These additional kaons come 
partially from double-0 production ifSOl with again small Q values since the four-K mass spec- 
trum has a steep threshold enhancement in the mass range from 2. 1 to 2.3 GeV [SOJ . This would 
mean that the effective contributions, Fig. ll22[ could increase by as much as a factor of 1.5, see 
below. 

Another candidate resonance for low-pT- kaon production is the A(1520) in the NK decay 
channel with its small Q value of 87 MeV. Measurements at ISR energy [j8l]| and at the SPS lf35l 
have been combined in Fig. ll23h to obtain an approximate dn/dxp distribution. 
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Figure 123: a) xp and b) distributions of A(1520) 

The full line in Fig. I123h is a hand-interpolation of these data that has been used in the 
Monte Carlo simulation. It integrates to the number (TiA(i520)) = 0.0219 per inelastic event or a 
total inclusive cross section of 0.697 mb. Since no data on the corresponding distribution are 
available the fiiio f^XF^Pr) = e~^^'^, with B = 2.9 as given in [[8T1 and shown in Fig. 1123b 
has been used. 

Due to the rather flat xp distribution of the A(1520) which is typical of neutral strange 
baryons, the ratio between decay and inclusive K~ shows a characteristic increase from xp = 
to a maximum at ~ 0.3 as shown in Fig. 11241 
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Figure 124: Ratio R^^^ = K^/Kj„^.j as a function of a) for different px and b) pr for different 

Xp 

Again, as for the decay, the very sharp enhancement of R^es towards low px is evident 
whereas the contribution to the inclusive kaon yield vanishes at about pr = 1 GeV/c. On the 
other hand the xp distribution of R^^ is complementary to the one from (p decay in its xp 
dependence such that the sum of the two contributions becomes rather xj? independent. This is 
evident in the combined ratio Rf^^ = K^_,_^/Kjj^^[ shown in Fig. ll25l where the yield from 
has been multiplied by a factor 1 .5 in order to make up for the production of additional kaons 
in (p production, see above. 
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Figure 125: Ratio Rf^^ = K^+A/^mci ^) ^ function of for different px and b) as a function 
of px averaged over the xp range < x^ < 0.4. In panel b) the error bars give the variation 
with Xp around the average. The full line represents the relative increase of the cross section 
as a function of y/s. Fig. U3h . normalized at = GeV/c 

As the combined contributions from and A(1520) decay reach about 25% of the total 
K~ yield at low px and y/s =17.2 GeV this discussion shows again the importance of resonance 
decay for the understanding of inclusive hadron production, in this particular case for the low 
Px behaviour of the kaon cross sections. This is also evident in the px distribution shown in 
Fig. 1125b which is very similar to the low px enhancement with ^/s shown in Fig. |73l The 
s-dependence in the region below p-r ~ 1 GeV/c will be determined by the s-dependence of 
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and Y* production with respect to the other contributing resonances. If the latter contributions 
rise as little as in the region 0.8 < < 1 GeV/c (p and Y* decays will become dominant in the 
ISR energy range. 

12.2 D(1860) decay 

In complement to the discussion of the low pr area of kaon production in the preceding 
section, it is interesting to look for resonance decay contributions in the high px region of 
Pt > i GeV/c. Here high mass mesonic resonances with sizeable decay branching fractions 
into 2 or 3 body final states including kaons will contribute. Although there is a large number 
of non-strange and strange resonances in the mass range above 1.5 GeV fulfilling this criterion, 
the charm mesons D°^(1860) will be regarded here as an example of heavy flavour production 
and decay. In fact the charm production threshold is crossed in the SPS energy range and the 
charm yields will start to saturate at p+p collider energies where beauty meson production will 
give access to still higher transverse momentum ranges. 

Close to 100% of all charmed meson decays end up in final state kaons, either in semi- 
leptonic or hadronic decay modes. Most of these are few body decays with large Q values, like 
Klz/,K*lz/ in the semi-leptonic and Kvr, Knn and Knnn in the hadronic case. Given the high 
D mass, the addition of one or two pions in the final state will not change the phase space 
distribution of the kaons appreciably. The two body decay of charmed mesons into K will 
therefore be studied in the following. 

One of the rare measurements of charm production in p+p interactions by the LEBC- 
EHS collaboration ll82l at the CERN SPS will be used to establish the input xp and px distri- 
butions as shown in Fig. 11261 
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Figure 126: a) dn/dxp as a function of xp and b) da /dp^ as a function of p^ of D( 1860) 

The inter/extrapolated dn/dxp distribution and the fitted Gaussian da /dp\ = Ae^^^-^^^T 
are shown in panels a) and b) as full lines. The integration of this parametrization yields 
{nj)0±) = 0.000944 per inelastic event corresponding to a cross section of 29 fih. This cross 
section, at y/s = 27 GeV, contradicts an upper limit of less than 10 /ib established from the 
study of muon pair production [|83ll at this energy. This discrepancy notwithstanding, the ef- 
fective cross section has been reduced to 20 /ib taking into account the steep s-dependence 
for the following comparison to inclusive K~ cross sections at the energy of the NA49 exper- 
iment, y/s = 17.2 GeV. The effective branching fraction of charm meson pairs into may 
be estimated from [|82l to about 47%. In order to take into account the softening of the decay 
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kaon spectra in multibody decays, a conservative value of 30% has been used for the following 
two-body decay simulation. 

Typical px distributions of the invariant K~ cross section from charm meson decay are 
shown in Fig. ll27l for two values of xp- These distributions are compared to the total inclusive 
K~ yields normalized to the decay distribution dXpr = GeV/c. 




[GeV/c] [GeV/c] p^ [GeV/c] 

Figure 127: Invariant cross sections from charm meson decay a) for xp = 0, b) for xp = 0.3 
as functions of px at y/s = 17.2 GeV. The corresponding total inclusive yields (dashed lines) 
are shown for comparison normalized to the decay distributions atpT = GeV/c. Panel c) shows 
the absolute ratio -Rcharm in percent as a function of px for the two xp values 

The decay kaons evidently show a much wider pp distributions than the inclusive K~. 
The relative increase of the ratio 

-Rchaim = „^ (38) 
Nncl 

is shown in Fig. 1127b as a function of pp for the two xp values of and 0.3. -Rcharm increases 
steeply with pp from values of less than 0.1% at low pp to more than \% dX pp = 2 GeV/c, 
whereas the ratio of the total inclusive cross sections is of order 0.15%. This increase will 
clearly continue at pp > 2 GeV/c. The situation is quantified for the complete xp/pp plane in 
Fig. [128] which shows -Rchami as a function of xp for different values of pp. 




Figure 128: -Rchaim as a function of for fixed values of pp, in % 
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A characteristic pattern emerges. At low xp -Rcharm gains about one order of magnitude 
between pt = and Pt = 2 GeV/c. This gain increases with x p and reaches more than two orders 
of magnitude at xp = 0.5 (see also the discussion of two body decays in [134|| ). In view of the 
sizeable experimental uncertainties still involved with charm production in p+p interactions, the 
percentage scale of the observed pattern should be taken as an indication rather than a precise 
prediction. Scale variations of up to a factor of two are easily possible should more precise 
measurements become available. It is the relative evolution with x p and pp which is unavoidably 
involved with heavy flavour decay given the precisely measured, large branching fractions into 
few body decays. Taking into account the rapid increase of the total charm cross section with 
-y/i there is no doubt that heavy flavour decay will represent an important contribution to the 
total kaon yields at large pp and at large xp already in the ISR energy range. 

12.3 Non-thermal behaviour of the decay products 

Transverse mass distributions of the inclusively produced kaons have been presented in 
Sect. 16.41 above. The inverse slopes of both K"^ and K~ show a strong variation with {nip — mK) 
from about 150 MeV at low nip — niK to 200 MeV at the upper limit of nip — niK available in 
this experiment. In this context it is interesting to have a look at the inverse slope parameters of 
the decay kaons from the 0(1020), A(1520) and D(1860) discussed above and shown in Fig. 1 1291 
as a function of nip — m^- 
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Figure 129: Inverse slopes of from the decay of 0(1020), A(1520) and D(1860) as a function 
of mp — niK- The result for inclusive K~ production is also shown 

Evidently the inverse slopes of K~ from 0(1020) and A(1520) decay cluster around the 
low "temperature" values of 80-100 MeV, whereas from charm decay show inverse slopes 
between 200 and 300 MeV. This is of course a result of the largely different Q values of the 
respective decays convoluted with the sizeable transverse momentum of the parent particles 
which gives them a mean transverse velocity {I3p) ~ 0.3-0.4. 

In thermal models such anomalies are not a priori foreseen, as all secondary hadrons are 
supposed to have Boltzmann-type distributions in mp with a unique inverse slope characteristic 
of the hadronic reaction involved. In Hagedorn's thermodynamic model for instance this "black 
body" radiation of hadrons happens from "fireballs" which are not allowed to have transverse 
momentum. Hagedom [[84l has in fact realized that decay products are non-thermal in the above 
sense for the decays A(1230) — )■ Nvr and p(770) — )■ tttt, albeit in a non-complete argumentation 
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as these parent resonances were in fact allowed transverse momentum and were taken as 6- 
functions at their PDG mass values (see 041 concerning the importance of the proper resonance 
mass distribution). If, as argued above, the majority of final state hadrons stem from the few- 
body decay of resonances which have important transverse degrees of freedom, the concept 
of a unique hadronic "temperature" in p+p interactions becomes an artefact. The fact that this 
concept is not able to explain the evolution of particle yields towards high transverse momentum 
and high ^/s, and, by the way, towards nuclear interactions without the introduction of ad-hoc 
concepts like the Quark-Gluon Plasma, see [|84l . has its origin in the same deficiency. 

13 Data summary 

After the detailed discussion of charged and neutral kaon yields in the preceding sections 
it is now mandatory to summarize the obtained results and to compare them to existing studies 
of global kaon production. The single differential, px integrated invariant cross sections F{xp = 
0), see Eq. [151 and the total yields elaborated in Sects. [10] and [U] are listed in Table[II]for K+, 
K-and K%. 
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Table 11: Single differential, px integrated cross sections F{xp = 0) in mb and total yields for 
K+, and for 31 values of y/s. The values in brackets for are derived from the cross 
sections and yields for the charged kaons under the assumption (n^o ) = 0.5{{n^+) + (^k-)) 



95 



A look at this Table shows that the Kg yields present by far the most dense and consistent 
coverage of the y/s scale from threshold up to collider energies, as compared to the results for 
charged kaons. This has already been evoked in Sect.[TTl see Fig. 1 1201 

13.1 Total kaon yields 

In a first attempt at establishing a consistent s-dependence from these data, the total 
yields {n^+), (n^-) and (n^o ) will be treated. These quantities are shown in Fig. 11301 
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Figure 130: Total yields {n^+), (t^k-) ^^id (^ro ) a function of y/s. The full line through the 
K° results is an eyeball fit, the lines through the K+ and data are derived from Fig. 11311 
below. The full circles in the data correspond to 0.5{{n^+) + {n^^-)) established at the 
corresponding ^/s values 

Whereas the coverage in the scale is dense and continuous for (n^o^), the corre- 
sponding data for the charged kaons show wide gaps in the range 7 < ^/s < 17 GeV and above 
^/s = 63 GeV. In this upper energy range, the extrapolation from ISR to RHIC energy has been 
evaluated in Sect. 110.41 The situation towards lower energies is confounded by the fact that the 
available data at ^/s = 4.9 and 6.8 GeV are evidently doubtful by internal inconsistency, see 
Sect. 110.21 The following procedure has therefore been followed to come to a consistent de- 
scription of the s-dependence. In a first step an eyeball fit through the K^. data is established 
over the full y/s scale, see the full line in Fig. 11301 This fit gives a consistent description of 
the situation within point-by-point variations of typically 10-20%. In a second step the ratios 
(nK+)/(^KO ) and {nj^-)/{nT^a) are obtained from the available data, see Table [TTl These ratios 
are presented in Fig. I131[ 
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Figure 131: Ratios (^k+)/(^k?,) (^k-)/(^k|^) ^ function of ^/s. The full lines are eyeball 
fits through the data at < 4.9 GeV and > 6.8 GeV 



A smooth y/s dependence is imposed on these data points between the low energy range 
y/s < 4.8 GeV and the higher energies y/s > 6.8 GeV/c using the fact that the cross sections 
in the PS energy range have been shown to deviate upwards (Sect. 110.21) . The (?t^k+)/(^ko ) 
(nK-)/(nKO ) ratios thus obtained are then used to produce the smooth lines in Fig. 11301 through 
the K"*" and the data by multiplying with the Kg interpolation. 

It is interesting to compare these interpolated results with the global study of kaon yields 
by Rossi et al. lIH which dates from 1975 but is still widely used today [[85l . Their fits to the K"*" 
and data are shown in Fig. ll32l in comparison with the present results. 
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Figure 132: {n^+) and {n^-) from ref. as a function of y/s, full lines, compared to the 
interpolated K+, K~ and yields Fig. 11301 (broken lines) 
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There are evidently major deviations over the full y/s range. Especially flagrant is the 
fact that the yields of [|4l are above the multiplicities between y/s =10 and 50 GeV. This 
is clearly unphysical. The relative deviations between the two attempts are shown in Fig. 11331 
on a percent scale. 
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Figure 133: Relative deviations between the present study of charged kaon yields and JH in 
percent as a function of y/s 



These deviations are in the range from -40 to -1-80% and demonstrate which order of 
magnitude of systematic effects must be expected if using existing parametrizations. The K 
analysis of [jUl is particularly questionable in the range 8 < ^/s < 12 GeV by using proton- 
nucleus data from a Serpukhov experiment If86]| with protons on Aluminium in an x^? range 
above 0.3 and for only 3 fixed lab angles between and 12 mrad. It remains a mystery, con- 
sidering the poor state of knowledge about proton-nucleus interactions and their normalization, 
especially in the strange sector, how these data could be translated into total kaon yields in p-i-p 
collisions. 



13.2 Pt integrated invariant cross sections at a; jr = 

As far as the px integrated invariant cross sections at xp = 0, F(xf = 0) are concerned, 
the experimental situation is similar to the one for the total kaon yields with the exception 
that only very scarce data below y/s = 6.8 GeV for charged kaons and below 4.9 GeV for K° 
are available. Again the K^. data may be used as a reference in establishing a consistent y/s 
dependence as shown in Fig. 11341 

The full line through the K^. data is an eyeball fit down to y/s = 4.9 GeV. As the UA5 
data [69] and the ISR extrapolation at y/s = 200 GeV show a 26% difference, see Sect. 110. 6?T1 
the fit has been chosen to pass 13% above the data extrapolation and 13% below the UA5 data 
which allows a smooth continuation to y/s = 540 GeV. In order to obtain a fit through the 
charged kaon data, a reference to K'g has been used by plotting the ratio between the charged 
kaon and interpolated data as presented in Fig. 11351 

These ratios may be smoothly connected by an eyeball fit between y/s = 200 GeV (both 
K"*" and shifted upwards by 13%) and y/s =11.5 GeV. As already visible for the case of 
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Figure 134: F{xf = 0, ^/s) as a function of a/s for K° (open circles: direct measurements, full 
circles: obtained as average from K+ and K~), K"*" (squares) and (triangles) 
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Figure 135: Ratio between K+ and interpolated data (full circles) and between K and 
interpolated data (open circles) as a function of y/s 



total yields, Fig. 11311 the data fl\ at ^/s = 6.8 GeV do not fall on the smooth extrapolation 
below ^/s = 1 1.5 GeV indicated in Fig. 11351 with a deviation of about 60% for K+ and 25% for 
K . This complies with the discussion of these data in Sect. 110.21 In view of this the following 
attempt to nevertheless obtain an approximate description of F{xf = 0) at low energy has been 



99 



followed. The fit to the data together with the and K^/K^ ratios has been used to 

obtain the charged kaon cross sections down to y/s = 4.9 GeV. The line through the K+ data 
has then been continued to the single measured K+ cross section at \/s = 2.9 GeV. The fit for 
has been discontinued at y/s = 4.9 GeV. This admittedly rather daring procedure produces 
nevertheless a consistent overall picture with systematic errors below y/s ~ 5 GeV on the level 
of about 20%. Only new precision measurements in this energy region may help to improve on 
this unsatisfactory situation. 

13.3 K+/K- ratios 

The data interpolation described above allows also an overview of the K+/K^ ratios both 
for the total yields and for the px integrated cross sections at Xi? = 0. This is shown in Fig. 11361 
where the full lines refer to the eyeball fits in Figs . 1 1 30] and 1 1 341 and the data points to Table [TTl 
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Figure 136: -Rk+k- ratios for total yields (open circles) and F(xf = 0) (closed circles) as a 

function of y/s. The error bar at y/s = 200 GeV marks the range of (0) / F^- (0) of published 
results by the RHIC experiments 

One feature of these dependences is the rather slow approach to K+/K = 1 for 
-^K+(0)/-^K-(0) with increasing y/s. For y/s = 200 GeV the ratio of the pr integrated cross 
sections at x^? = is 1.14 for the extrapolated ISR data as compared to a range of values be- 
tween 1 .03 and 1 .08 published by the RHIC experiments Il25l - [30l . The reason for this further 
discrepancy might be again the low trigger cross sections at RHIC with respect to the total 
inelastic cross section which enhances central collisions and hence kaon ratios closer to unity. 

14 Conclusions 

The new data on the inclusive production of charged kaons in p+p interactions at SPS 
energy presented here complete a detailed study on charged secondary hadrons including pions 
and baryons in the framework of the NA49 experiment at the CERN SPS. These data offer the 
possibility to check the sum rules of mean charged multiplicity and charge conservation. It is 
demonstrated that both constraints are fulfilled within the tight error limits of about 2% as they 
have been quoted for the systematic uncertainties in the independent evaluation of inclusive 
cross sections for the different particle types. 
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The extended coverage of the data in the xplpx plane, from = to 0.5 and from 
Pt = 0.05 to 1.7 GeV/c, allows a precise study of particle ratios from the same experiment, 
both for K/tt and K/baryon ratios. A detailed comparison to existing data in the SPS/Fermilab 
energy range shows in general good agreement in the limited phase space regions available, in 
particular also for the complete set of particle ratios, with some exceptions essentially due to 
normalization problems. 

As the interaction energy of ^/s = 17.2 GeV is located at a strategical point be- 
tween threshold-dominated and scaling phenomena at lower and higher energies, respec- 
tively, a new and complete study of the s-dependence of kaon production including K" has 
been attempted using the new NA49 data as a reference. This study covers the energy range 
3 < < 1800 GeV and aims at establishing an internally consistent picture of kaon produc- 
tion as far as this is possible with the often restricted and contradictory available data. Through- 
out, the use of data parametrization with simple arithmetic formulations has been avoided in 
order to take the rather complex dependence of the measured particle yields on the kinematic 
variables fully into account. 

This study reveals basic weaknesses in the existing data base both at lower and higher 
V^. At PS/AGS energies the charged kaon data suffer from large systematic inconsistencies, 
and the almost complete absence of differential data at Serpukhov energies renders the estab- 
lishment of integrated yields hazardous to say the least. The extension into the ISR energy 
range on the other hand, using all available data, gives new insights into the complex evolution 
of strangeness yields as functions ofxp and, in particular, pr resulting in explicitly non-thermal 
transverse distributions. These findings are discussed in connection with some typical examples 
of resonance production and decay which are relevant to this phenomenology. 

In addition to the charged kaon cross sections, it has been found useful and necessary to 
also look at the evolution of the Kg yields. In fact the relation = 0.5(K"'" -|- K~) which is 
found valid within the experimental precision at least in the energy range ^/s > 5 GeV provides 
a strong constraint on the overall data consistency. In this context the early bubble chamber 
work up to ^/s 28 GeV proves to be essential due to its internal consistency and its superior 
precision concerning the overall normalization. 

As far as the extension of the study to the RHIC and p-i-p collider energies is concerned, 
a rather disturbing overall picture emerges. Evidently, the published results do not represent a 
decisive improvement as far as precision and internal consistency are concerned in comparison 
to the lower energy data which in most cases date back by more than 30 years. There are several 
reasons for this situation: 

- The study of soft hadronic production in elementary collisions is certainly not at the 
heart of the experimental programs at collider energies. On the contrary it is the discov- 
ery potential for "new" physics either in Heavy Ion interactions (RHIC) or within and 
beyond the Standard Model (p+p colliders) which defines the main priorities. 

- Precision studies of elementary hadronic production call for specific constraints both 
concerning accelerator layout and operation (as for instance vertex distributions and 
stability) and the experimental set-ups (trigger efficiency, data normalization, material 
budget at low total momenta). 

- The large multi-purpose detectors generally set up at the colliders are not really opti- 
mized for these constraints and small- size, dedicated experiments as they have been used 
in practically all the preceding lower energy work, are neither available nor planned. 

In view of these problems in a first step the extrapolation of the ISR results to ^/s = 200 
GeV has been attempted in order to obtain a common point of comparison between ISR ex- 
tension, RHIC and the lowest available p-i-p collider energy. At this energy the UA5 streamer 
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chamber data turn out to offer a reliable reference although the overall statistical errors are size- 
able. This is reminiscent of the bubble chamber work at lower energies which definitely benefits 
from the application of optical methods in terms of reconstruction efficiency and normalization. 
In contrast all RHIC data show large systematic offsets and a general weakness towards the 
lower cut-off in px which lies in general in the region 0.4-0.7 GeV/c. 

This new study of the s-dependence of charged and neutral kaon production results in 
smoothed interpolations of the central, px integrated invariant yields F{xf = 0, a/s) and of the 
total kaon multiplicities as they are presented in Sect. [13] of this paper. 
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